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Abstract 

The non-stationary Kinetic theory of Israel and Stewart 
[16] is extended from the case of a gas consisting of only 
one species of particle to the case where the gas consists 
of an arbitrary number of species. A particular species may 
have non-zero or zero rest mass (e.g. electrons or photons). 
It will be shown that generalized fitting conditions, which 
determine the temperature and relativistic chemical 
potential of each species, are required to adequately treat 
the physics of the gas, which nevertheless will be 
independent of the choice of the fitting conditions to first 
order. Equations which describe the transport of thermal and 
viscous effects will be derived; and coefficients of thermal 
coductivity, bulk viscosity, and shear viscosity will be 
defined. We shall also obtain the entropy production for any 


particle species. 
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De i minoauction 

Irreversible relativistic thermodynamics has undergone, 
in the last fourteen years, an extension of its 
applicability. Originally, it was a stationary theory in 
which gradients of the macroscopic deviations from 
equilibrium, such as heat flux and viscous stresses, were 
considered negligible on the scale of the mean free path. 
However, this theory had one great defect. It predicted 
infinite speeds of propagation of thermal and viscous 
effects whereas we would expect such speeds to be about the 
mean molecular speed and certainly not greater than the 
speed of light (Israel[14]). 

Early attempts at resolving this problem were focused 
on the heat transport equation (Fourier’s law) and consisted 
of the addition of ad hoc terms to convert it from a 
parabolic to a hyperbolic differential equation (infinite 
versus finite propagation speeds respectively). The 
suggestions of Cattaneo [5] and Vernotte [39] were further 
discussed by Kranys [20]. These authors proposed that a term 
proportional to the time derivative of the heat flux be 
added to Fourier’s law and showed that the resulting 
equation was hyperbolic with finite propagation speeds. 

In a paper on non-relativistic thermodymamics, Muller 
[28] suggested that the expression for the entropy was 
incomplete and, in addition to the number density, pressure, 
and energy density, must include as independent variables 


the heat flux and viscous stresses. This suggestion is in 
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fact the foundation of the successful relativistic theory 
and is essentially an extension of the theory from the 
stationary to the non-stationary regime. Muller went on to 
show that hyperbolic equations of transport resulted as a 
consequence of this formalism. 

Subsequent discussions of the propagation speed problem 
focussed on the non-stationary aspect of the relativistic 
theory via the relativistic Grad method of moments. Kranys 
P21, 221 extended Chernikovss 17-8) thirteen moment 
stationary theory to the non-stationary case and obtained 
hyperbolic propagation speeds. Stewart [33] showed that by 
retention of non-stationary terms in the fourteen-moment 
theory, the upper bound on propagation speeds was 73/5 c. 
Discussions of the fourteen moment case were also presented 
independently by Marle [27] and Kranys [23]. 

The discussions cited above were in the realm of 
Kinetic theory. Within the context of relativistic 
phenomenological theory, however, the problem of a causal 
theory which predicted finite propagation speeds was not 
resolved until Israel [14] independently rediscovered and 
extended the idea of Muller to the relativistic case. This 
theory was later combined and compared with the fourteen 
moment Kinetic theory by Israel and Stewart [15,16]. An 
analysis of the propagation speeds of thermal and viscous 
effects by an examination of the characteristics of the 
transport equations was reported by Stewart [34] and Israel 


and Stewart [17]. An analysis via Fourier analysis of the 
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transport equations was performed by Kranys [24]. The upper 
pound. on the Speed of propagation, of thermal and viscous 
effects appears to be established as 73/5 « [16]. The 
extension of this theory to polarizable media in the 
presence of electric and magnetic fields has been performed. 
We shall not discuss that theory herein but refer the reader 
to a recent article by Israel and Stewart [18] which 
performs this extension and also discusses the background of 
this theory. 

In summary, we have at the present time a 
phenomenological theory for relativistic simple media and 
mixtures which takes into account transient effects and 
predicts finite propagation speeds of thermal and viscous 
effects; we also have a relativistic transient Kinetic 
theory for single species gases. There are astrophysical 
situations, however, where this latter theory is not 
applicable because we have to deal with plasmas. For 
example, such situations are: the accretion of matter 
through the boundary between a star and a hypothetical 
neutron star or black hole in the star’s core [38]; 
accretion disks around neutron stars or black holes; and the 
leptonic era in the early history of the universe. In these 
situations we imperatively require an extension of transient 


performed in this thesis. Furthermore; we apply the extended 


theory to a case of special astrophysical interest: a 


mixture of matter and zero-mass particles (photons or 
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neutrinos); that is, we study relativistic radiative 
transport with the inclusion of transient effects. 

To illustrate the need for transient theory as opposed 
to quasi-stationary theory in astrophysical situations, we 
shall briefly discuss accretion disks and the leptonic era 
in the early history of the universe. Consider a black hole 
which is accreting matter. If the accreting matter has a 
large amount of angular momentum then an accretion disk is 
formed [4]. The principal model discussed in the literature 
is the thin disk model where the half thickness of the disk, 
h, is much smaller than the distance to the black hole, r, 
i.e. h/r << 1. The matter in the disk moves in Keplerian 
orbits and viscosity between adjacent rings transports 
angular momentum outward and matter inward. The viscosity 
also acts as the principal source of heating in the disk. 
Most of the radiation is produced in a small region near the 
black hole and ionizes the infalling matter to distances 
beyond the accretion radius of the black hole. Thus the disk 
consists wholely of a highly ionized plasma. The disk has 
three zones of interest: the outer zone where the pressure 
is dominated by the matter pressure and the major source of 
opacity is due to free-free collisions (Bremsstrahlung); the 
middle region where gas pressure still dominates but the 
major source of opacity is due to electron scattering; and 
the inner zone where the radiation pressure dominates and 


opacity is due mainly to electron scattering. 
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The vertical structure of the disk is determined by a 
balance between compressional tidal forces and the outward 
pressure and radiation flux. The optical depth t is given by 
t= h/XA [25] where A is the mean free path of the photons. 
iheshalt -thichnessvof iherdiska Inn @Se audiStance sca levor 
instabilities: ands? luctuations im the disk [25,321). Normally 
a stationary radiative diffusion equation is employed to 
help find the vertical structure of the disk. In the 
Opts calilVerthick case. a>) (and: anethe optically thin case, 
(RE this equation is valid; however, when t=1 that is, the 
disk is neither optically thick or thin, then we are in the 
non-stationary regime and we must use transient theory. 

Novikov and Thorne [30] have developed a model of a 
thin disk. In this model they note that in the outer and 
middle zones, the disk is optically thick. However they 
specify a solution in which a small region in the inner zone 
is optically thin. In this model therefore, there must be a 
transition zone between the optically thick and optically 
thin regions where we must employ transient theory to 
determine the vertical structure. Eardly et al. [9] have 
also constructed a disk model in which they find the inner 
region of the disk to be marginally optically thin to 
electron scattering. Hence we should apply transient theory 
to this model as well. Other models may also require 
transient theory because thin disk models are secularly 
unstable in the inner zone so that the thin disk model may 


be invalid in the inner zone [26]. The disk problably 
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becomes, in this region, a small cloud around the black 
hole; In this case the optical cepth could be close to, one. 
Also, convective turbulence may cause the disk to thicken so 
that the half thickness of the disk could become the same 
size as the photon mean free path [19]. We conclude 
therefore that transient theory is a valid or a necessary 
approach to accretion disk physics. 

The history of the universe is usually divided into 
four eras: the hadronic era when the temperature, T, was 
greater than 10!* °K; the leptonic era when 
hot ek > 1 Silo?’ ek) the radiation dominated era when 
103° eke ies el0 SK and sthe matter dominated era when 
Te =10) | oR elhise last tera comprises 99994 of the history of 
the universe. We restrict our attention to the leptonic era. 
At this time the universe consisted of a mixture of 
electrons, muons, neutrinos, photons, and their 
anti-particles. We can estimate the temperature when 
transient effects become important by comparing the rate of 
expansion of the universe, H, to the rate of interactions , 
a where Ok is the weak interaction cross-section and 
n is the lepton density. As the universe expands the 
particles will tend to go out of equilibrium with each other 
but their interactions tend to restore that equilibrium. 
Equilibrium will be maintained when Oj n/H >1 (rate of 
expansion is less than rate of interactions) but equilibrium 
will not be maintained when O 4 n/H < 1. We conclude that 


transient theory will become important when o ,n/H = AA, 
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Weinberg [43] has computed this ratio: 


Owk™ -{ T } exp{ - 20% x } ean 
‘ 1010 °x : . 


This ratio is unity for T = 1.3 x 10!! °K, Hence transient 
effects will become important when the temperature of the 
universe drops below this value. 

We shall now summarize the phenomenological theory for 
a simple (single component) fluid in a gravitational field; 
this summary allows us to illustrate the principal features 
of the non-stationary theory and to identify the essential 
difference between quasi-stationary and non-stationary 
theory. An arbitrary state of the fluid is described by 


three primary variables. These are the entropy flux 5°, the 


energy-momentum tensor poe and the number flux N° . The 
number flux and the energy-momentum tensor are assumed to be 
conserved and the entropy production is assumed to be 


positive: 


where the single stroke denotes covariant differentiation. 
The equilibrium state of the fluid (denoted by a 
superscript 0) is characterized by four properties. First, 
the entropy production is zero: S41 0. Secondly, there 
exists a unique time-like unit four vector u” which we call 


the flow vector, such that we have 


that is, the entropy flux, number flux, and the 
energy-momentum tensor are spatially isotropic. Here S is 
the entropy, n is the number density, p is the energy 
density, P is the pressure, and A eines eee ys the spatial 
projection operator of uv Thirdly, each equilibrium state is 
eharacteri zed by an equatiom @f state s=sm@m.p) which 


determines the entropy and from which we can find the 


pressure by 
S= (9 +°P)/T ran , Gi) 


where T is the temperature anda is the relativistic 
chemical potential of the equilibrium state. Finally, the 
flow vector u- is shear free and expansionless; and the 
relativistic chemical potential is constant: 
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The covariant formulation of the equilibrium form of 


the entropy flux may be expressed as 


a 


‘on 


where Wenn Equation (1.6) implies that 
a au ° 
= - adN = yal . Ld en) 


where the differentials are constrained to displacements 
between equilibrium states. 

To investigate non-equilibrium states we release the 
constraints on the differentials in equation (1.7) and 
assume that they apply to displacements from equilibrium 
states to arbitrary nearby states: 


qe 2p. Gene = pvate : 


Here a and 6. are the variables of a nearby equilibrium 
state such that the deviations from equilibrium eee and 
eb _ pk are small quantities compared to N° and 8 and 
are said to be of first order. Equation (1.8) suggests that 


equation (1.6) be generalized to 
steep eat = oN: oem te Ge (479) 


where P is the pressure of a nearby equilibrium state and 
gh is some unspecified second order term. 

Now the essential difference between quasi-stationary 
theory and non-stationary theory resides in the treatment of 
o- = In general, o will depend in some fashion on the 


macroscopic deviations from equilibrium. If these deviations 
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from equilibrium vary over large distance scales compared to 
the mean free path, then their space-time gradients are 
negligible to first order and hence ql, will be negligible 
to second order. Therefore aq, will not contribute to the 
entropy production which is calculated to second order; and 
we may therefore neglect QF an equa tion.( two), a 
non-stationary theory, however, the deviations from 
equilibrium vary over a length scale comparable with the 
mean free path. Consequently their space-time gradients are 
not negligible to first order and then Qi, is not negligible 
to second order. In this case, Qh, contributes to the 
entropy production and neglect of ae ie equation (1.9) is 
not justified; that is, in the non-stationary theory, Qh 
must be retained. 

In equilibrium, the unit time-like vector parallel to 
the number Wh, ce, and the unit time-like eigenvector of 


the energy-momentum tensor, coincide. In non-equilibrium 


Qa 
Une 
they differ by a small angle of first order. We can choose 


an arbitrary velocity u- within a small cone whose angle is 


of first order and includes Be and uy . Then the number 
density, the energy density, the pressure, and the entropy 
are independent of the choice of u- to first order. 
Furthermore, for a given choice of u' we can decompose the 


number flux and energy-momentum tensor uniquely: 


NY = nul + 44 , ju, =0; 
me eS itt ee TV AveeEeh ike ebeee qo” , (ae dO) 
uh” ="01 5 a =50.., Ty or ()) =: ne wu” + gh’ 


Here ;° TS Oeieor Culescirgtaes h is the momentum flux, 7 is the 


bulk stress, me are the viscous stresses, and Ace 


is the 
spatial projection operator of ihe pDartiiele drift and 
momentum flux change in first order if the choice of wu” 
changes. Consider however, the heat flux q. which we define 


as the energy flux relative to the particle flow; it is 


given by 


Then if the choice of u’ changes, q is unchanged to first 
order. 

A crucial aspect of transient phenomenological theory 
1s hes spect ication oF a in equation (1.9). To obtain 
linear phenomenological laws it is sufficient to assume that 


Y is a quadratic function of the macroscopic deviations 


@) 
from equilibrium. 7, q etc. In particular, the simplest 
approximation which leads to hyperbolic phenomenological 
laws is the hydrodynamical description which assumes that an 


arbitrary state of the gas close to equilibrium can be 


specified completely by the number flux and the 
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energy-momentum tensor. Hence, the entropy flux is a 
function solely of the number flux and the energy-momentum 


tensor. In this case q is given by 


al B 


HL u D r a 
ali = =< 
OD RSS ee eaeiysl Se a a) (4.12) 
. ee A Wel ale ak A U A Ww 
nm a,49 7, + Parca ( 5 h hu +h 7, oe 
where Cy» Oy, Bo» Bo> and B., are undetermined 


thermodynamical functions. We require that HQ /s0 which 
implies, for all states characterized by a given number and 
energy density, that equilibrium has the largest entropy; 
furthermore, it can be show to imply that the fluid is 
dissipative, that is, it has positive relaxation times. 

The phenomenological laws (transport equations) may be 
inferred from the requirement that the entropy production is 
postive, S] 420 if we assume there -is a linear relationship 
between 1, qj: and 1°® and their gradients, the gradients of 
the thermodynamical variables and the flow vector, the 
shear, and the volume expansion. Consequently with the 


choice =u, we obtain 


o : eg 

US ( “E| % Bot “4 lee : 

qh = KTAuA( nt Q B q =O ee ap Chk Wi y 5 : (le ie, 
perp (2 ey Ouale eel yee ae 

Une = 27. ( rculys + BoM Wy - Cad oels ee: 


where by is the bulk viscosity, x is the thermal 
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Gonduelivirty ; So 18 shear viscosity, and the angular brackets 
denote the trace free spatial part of the tensor enclosed by 
them. If we had chosen w =u, we would have obtained a set of 
equations similar to those above but with some of the 
coefficients slightly changed. These equations are 
hyperbolic and predict finite propagation speeds. 

The Kinetic theory of a single species gas duplicates 
the features presented above for the phenomenological 
theory. Briefly, Kinetic theory postulates a distribution 
function for the gas and the Boltzmann equation which 
governs the distribution function. Then, in terms of this 
distribution function, the number flux, the energy-momentum 
tensor, and the entropy flux can be defined. The requirement 
So for equilibrium uniquely specifies the distribution 
function as a function of the temperature, the relativistic 
chemical potential, and the flow vector u' . The transition 
to non-equilibrium states of the gas close to equilibrium is 
achieved via the relativistic Grad fourteen moment 
approximation method, which assumes that the distribution 
function is expressable in terms of an equilibrium 
Grstributmion wtunetion and atins @ order.quddratic function 
of the particle four-momentum. We can then find the 
distribution function ine terms of the: structures oiven vin 
equations (1.10) and (1.11). Consequently the entropy flux 
is found to have the structure given by equations (1.9) and 
(1.12). The transport laws may be derived and come from the 


conservation laws for number flux and energy-momentum and 
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the balance equation for the third tensor moment of the 
distribution function, which we herein call the 
double-momentum flux. An important aspect of the kinetic 
theory as opposed to the phenomenological theory is that we 
can actually specify the functional form of the transport 

Q 


coefficients and relaxation times a ane eat icl B, for 


0° 
a gas. 

The problem considered in this thesis is the extension 
of the relativistic kinetic theory, as presented in the 
formulation of Israel and Stewart [16], from the single 
species gas to a gas which contains many particle species. 
We do not wish to restrict the analysis to species which 
have only non-zero rest mass, but also wish to consider 
species which have zero rest mass, such as photons and 
neutrinos. The theory of Israel and Stewart is suited only 
for non-zero rest mass particle species; however, we shal] 
extend the theory in an appropriate fashion to handle zero 
mass particle species. 

One might expect that we shall find the structure of 
the entropy flux to be similar to equations (1.9) and (1.12) 
for each component of the gas. Futhermore, one might expect 
that we should obtain transport laws that resemble equation 
(1.13) although of a more complex nature. Also, one might 
expect that there are more functions to be specified like 
Qa, Gas etc. Indeed, the derivation of all of these details 
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We shall briefly outline the format of the thesis and 
the major points of each chapter. The main body of our 
thesis consists of six chapters plus two appendices. Chapter 
Il consists of the fundamentals of Boltzmann Kinetic theory. 
This consists of specifying the Boltzmann equation and 
deriving the master balance equation, special cases of which 
are the mass and number fluxes, the energy-momentum tensor 
and the entropy flux. Chapter III applies this theory to the 
equilibrium situation. The two important aspects of this 
chapter are the introduction of the thermodynamic functions 
necessary for all of the subsequent analysis of the physics, 
and the equilibrium structures of the number flux, the 
energy-momentum tensor, and the entropy flux. The material 
presented in these two chapters is not new but its 
presentation is necessary to provide the foundation upon 
which the subsequent analysis relies. 

In chapter IV we analyse the non-equilibrium situation 
of the gas via the Grad fourteen moment approximation. It is 
here that we shall introduce the idea of fitting and frame 
changes under which the mathematical description of the 
physics will be shown to be invariant to first order. Also 
the deviations from equilibrium are solved for, in terms of 
physical quantities such as heat flux and the viscous 
stresses, by a method which differs from the approach 
normally used (Israel[12], Israel and Stewart [15,16]). This 
approach leads to a more aesthetic and compact notation. In 


terms of this solution, we shall then obtain an expression 
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for the entropy flux. 

Chapter V examines the non-stationary (transient) 
aspects of the theory. It is here that we compute the 
derivatives of physical quantities such as the mass flux and 
the energy-momentum tensor. We also compute the derivative 
of the entropy flux for any particle species, that is, we 
compute the entropy production for each species in the gas. 
These computations proceed by a different, although 
equivalent, approach than that of Israel and Stewart [16] 
for the single species case. Furthermore, the equations of 
transport of thermal and viscous effects for non-zero rest 
mass particle species are derived; and the coefficients of 
thermal conductivity, buik viscosity, and shear viscosity 
for each species are defined. 

Chapter VI considers the problem of extending the Grad 
method of moments to zero rest mass particle species such as 
photons and neutrinos. The solution for the deviations from 
equilibrium in terms of physically meaningful quantities 
leads to, when inserted into the Boltzmann equation, the 
non-stationary transport equations. Finally, we examine the 
entropy flux and define the coefficients of heat 
Conductivity, shear vascosity,, ama bulk =Viscosily for zero 
rest mass particle species. The scenario envisioned in this 
chapter is a mixture of matter and radiation, that is, this 
chapter deals with cnaaeacet radiative transfer. We noted 
earlier that this situation is of special astrophysical 


interest and would be applicable to important situations, 
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e.g. accretion disks and the lepton era of the universe. 

Chapter VII deals with the computation of quantities 
which depend upon the collision cross sections and the 
deviations from equilibrium; these quantities are necessary 
for the solution of the transport equations. This chapter 
should really be regarded as a Super appendix because it 
deals with a very complicated analysis to provide minor 
although necessary information for the discussion of the 
DiVS1¢s. 

The first appendix deals with the functional form of 
the thermodymanic functions introduced in chapter III. 
Various relationships between these functions are cited and 
some other frequently used results are stated. 

The second appendix is a list of the symbols employed 
in this thesis. It is arranged alphabetically, Latin letters 
first followed by Greek letters. The list is necessary 
because so many symbols are used that the reader may lose 
track of the meaning of a particular symbol. Instead of 
trying to find the meaning in the text, all he or she has to 
do is refer to this list for a quick reminder. 

As a final note, we declare our space-time convention. 
In a local Lorentz frame the metric takes the following 
form: 

Mag ily = HEE Cia iby 2 (4.14) 
The space-time coordinates x” have the Euclidean-Minkowski 


form 
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As a consequence of this convention, four vectors are 
time-like, null, or space-like if their lengths are 
negative, zero, or positive, respectively. Finally, 
throughout this thesis, we shall adopt units such that the 


speed of light, c, and Planck’s constant, h, are unity 


(h=c= le 
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Il. Boltzmann Kinetic Theory 

All of the calculations in this thesis are performed 
within the context of Boltzmann kinetic theory. Before we 
can examine the main topics of this thesis we must present 
the notation employed and the fundamentals of the kinetic 
ChHeernyjetiltseoSmeune DURDOSe Gime Lil > Cacbters 

First we shall present the notation employed and some 
basic definitions. Then we shall define the distribution 
Function Ny and present the Boltzmann equation as the 
equation which governs the evolution of the distribution 
function; at this point we shall also discuss the role of 
collisions of the gas particles. Once we have specified the 
Boltzmann equation we may derive from it a master balance 
equation which is then applied to find the balance equations 
for suitably defined quantities such as energy-momentum and 
entropy. Finally we discuss how these quantities may be 


decomposed with respect to an arbitrary comoving observer. 


A. The Boltzmann Equation 

Let us consider a multi-component gas. Each component 
of this gas is a particle species which we will identify by 
an Arabic number. These numbers will be generically 
represented by capital letters, for example, A = 1,2,3 etc., 
and used as subscripts on our mathematical symbols. 

A particle which is a member of species A will have a 
mass m , a charge 


A 
distinguishing characteristics) which are common to all 


e,, and a spin s, (plus any other 
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Dariicles oF species A. Each particle has 1ts own world-Jine 
ee) where T, 1s the world line parameter; we assume that 

a is an affine parameter for geodesics. Also, each particle 
has its own velocity w(t) = dx / at, and Kinetic momentum oe 


°o 


Qa Qa (0! 
= = = = 
We shall denote Pa = mW, where mM =m, WaWag for non-zero 


rest mass particles and Byes le Shoe 0 for zero rest mass 
particles. We shall refer to non-zero rest mass particles as 
massive particles; also, we shall refer to zero rest mass 
particles as massless particles. 

At each space-time point x we construct a 
four-dimensional momentum space, whose co-ordinates are the 
components of momentum Bi, , and which is a subspace of the 
eight-dimensional phase space BGasry) for species A. Not all 
of this momentum space is accessible to a particle. The 
condition PaPyom “me confines the particle to the surface of 
a pseudo-sphere in the momentum space. We say that the 
Dartiele 1s on—shel |. 

We choose, at x” , an arbitrary space-like element of 
three volume d= with surface normal nn. The element of 
volume of the accessible momentum space dv, is the four 
dimensional volume element yaa d'D, multiplied by the 
on-shell condition S(pop,t m2) O(p4) (eis 

The invariant distribution function for species A, 
Wa, D5) is defined by stating that [11] the number of noni 
lines of particles of species A which cross ay in the 


positive sense of the unit normal with momenta in the range 


of dv, is given by 
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Hy Caan (nw) (n.n“)dZdv, ; boee 


We postulate that the behaviour of N,(',p,) , and 
hence of the gas, is described by an equation of continuity 


in phase space with sources and sinks: 


Here N, 18 the abbreviated form for the distribution 
function, || means covariant differentiation holding 
momentum fixed by parallel propagation, a/8p, means 
differentiation with respect to momentum holding position 
fixed, and bp, /6t is the rate of change of momentum along a 
particle’s world line. Equation (2.2) may be regarded in a 
different manner. The left hand side is just the Liouville 
operator acting on the distribution function, and describes 
the rate of change of N, along the streamlines in phase 
space [13]. In the absence of collisions this rate of change 
is zero. The right hand side is therefore a correction to 
the Liouville equation to account for collisions. Equation 
(Qe) iseoumuscit zmannieeduation, 

Unless otherwise stated, we shall assume that there are 
no external fields other than gravity, so that the rate of 


change of momentum between collisions is zero, that is 


ép,/dt = 0. 
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Consider two particles of species A and B which undergo 
a collision. The incoming particles have momenta oe and sg 
and the outgoing particles have momenta P and Pe . The 


iiae *a *Q 
collision conserves momentum: Dies 


20 xa, Byte 
STDs SP Be The transition 
probability for this collision is denoted by W(D, Pp lPgsPp)« 
ihis transition probability ts a scalar funtion of 


* * 
Die a aun shar and has the property that [40] 


* * * * 
WOE Per, = W(Pp>Py [Pp oPq) : eae) 


We shall also assume that the transition probability 
satisfies the bilateral normalization property which appears 


here in a form suggested by Weinberg [42]: 
S * A A ay au 
Cheeta A Adv a 
= | WEP, sPpIPy Pp) Anda, . 


In equation (2.4) we have Me ete, N With e -2 for zero 


rest mass particles, g,=2s,+1 for non-zero rest mass 
particles, and ene nels classical particles; e=1 for bosons, 
ell for fermions, and al Tore classical ’parlicles. or 
classical particles, equation (2.4) reduces to the form of 
the bilateral normalization as stated by de Groot et. al. 
[11] and can be derived from the unitarity of the scattering 
matrix in quantum mechanics [40]. We note that equation 


(2.4) is more general than the assumption of detailed 


balancing since the latter immediately implies the former. 
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We can now write out the source and sink terms for the 
Boltzmann equation. Unless otherwise stated, we shall only 
consider binary collisions, which automatically preserve the 
number of particles of each species. The source terms are 
the number of particles created in a momentum cell by 
collisions; and the sink terms are the number of particles 
destroyed by collisions in that cell. The difference between 
the two terms (source minus sink) is D ces 

* x x x x k oO* 
eolienes L { Ny Ny W(Py »PglPy>Pp)A,AnaVaV RV, 


a ko 
) | Ny N, Wp, Ppl Py Pp) A ARIV,AV, AV, 


The a terms in equation (2.5) account for the Fermi 


exclusion and stimulated emission quantum effects. 


B. The Conservation Equations 

A master balance equation can now be derived from 
equation (2.2). We let & = ¥,d®,/dN, where ¥, is any tensor 
function of position or momentum and o, 1s a function of Ny 
alone. We multiply equation (2.2) by E, and integrate over 
the accessible momentum space: 
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After some algebra, we obtain our master balance 


equation 


" 
A | A a : ‘i = 
et. 2} ates Ante ite 


ow bet els 


vires fh sete 


pl ev aL SLE’ 
1°90 SO VLCe ‘7! 
| |) 4a thy Mer: 

’ inset 

i | ‘ nt , ar 
j ! 

, a | 

* 4? 
15, OV \ 
arg ' a 6G 
fa 


Ke i 
i 4 Mu ODT 
Fy Ts hu vii ai 
t 
* , | } } 
oy 7 5 Pe 
Ye ia-rOr Twa) 4 2 


swe stasgeale oth 


phe 


Be heroiic: 28s 


non Marri 


5e57- = Te 
i oaks 
a fif Delhesao- eure a ib "4 a] 
{ | -<Neet tr on ied 04 tt =) 
' | j ? : rept a] \ ime ca Geena b 
— a 
al j 
a rE writ wu ow! et 
2 as 
a nt | » 
* : ! iy oa)'s : 4 7 \ ct P > : 
: d 


J 


ia 


Sanne oct 


Ty ea 


ed 
me ro" Nisbaciem re 
ein +1 eaa! ve ie; aie as aia 
lbh vata Wily te saan 


" 
a, —_ = 
ty war Phd wi os 


$ i Biba) Fikes fT, J 


Lam 6S 


f ‘ - i's } 
a . . 
‘ tai 3 be Sila cot law 
hin : whiny ; ree | (toe ie ci Ate) ' 
ie biel Te ition oor 


ny 
7 oe fe” Wilt ree 4 tof af 
a 


| > 5 eee 
> ji oe Fed 


24 


6Y 
a ae A : 
{| ov satay, b = | Saran Cla | Sn Done hile (FA) 


where the single stroke on the left hand side means 
covariant differentiation. 
The expanded form of the last term in equation (2.7) is 


given by 
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where d'tv = dv dv dv dv... 
When we relabel the variables in the first term in 


equation (2.8) we obtain the following expression: 
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From equation (2.9) we see that if EAH ey Petia t is, & 
is a collisional invariant, then the right hand side of 
equation (2.9) is zero. If > is not a collisional 


invariant, we sum over all species so that equation (2.7) 


becomes 
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Ime Vast term in equation (2.10) becomes, via 
relabelling of variables in equation (2.8) and using 


equa lOneueses. 


ye ees Nghia es ) [EJN NAY ec Be (26113) 
ASCODLe AS ok ” Nip aoe an x ; 
A AB 
ees : : : 3 

where el = Ete, ~ 5, S Lt 1S a Summational invariant of 
Les Coll si6M, then del .= "0: 

The familiar conservation laws are now just particular 
eases Of equations (2.7) dndsi2-10) We shal) discuss each 
of these in turn. 


The number f1ux 1 is defined by 
= | nye dv : (2.12) 


bgt OPEN and ie alee equations (2./))-and (229) ‘give us the 


A 
number flux conservation equation: 


Qa 

“iil . (213) 
We note that this result depends heavily upon our decision 
to consider only elastic binary collisions. 


The mass flux M. is defined by 


A 
mi = | 1,5 dv (2a14) - 
A AA AGKs 
We let o,=N, and Y= my i consequently equations (2.7) and 
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(2.9) give us the mass flux conservation equation: 
QO — 
eae 5 C2154 


The energy-momentum tensor of particle species A is 


defined by 
ONS a 8B : ; 
r, = | NAPAMa dv, ; 2 109 


and we define the total energy-momentum tensor by 


Noting that [p,] = 0 for collisions we select ®, = N, and 


Weeih _ Then equations (2.10) and (2.11) give us the 


energy-momentum conservation equation: 
=O | omen) 


In contrast to this, equation (2.7) gives us the following 


result: 
a8 e a, 
oe ge | By Daas hie (2.19) 
A useful physical tensor denoted by a , and which we 


shall call the double-momentum flux, is defined by 
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Un = | NaPawaWy dv, : (2a 0) 
With the choices ® = N, and ¥, = paws , equation (2.7) gives 


us the fol lowing result: 
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The entropy flux for species A is defined by 


Oe = 2 o ; 
S, = 7k [ Ov, 04/8.) e,A,in(A,/g,)+(ex-1)N, Iw av, a en 2s 


and the total entropy flux is defined by 
Qa Qa 
S =) s : eee) 
i A 


where kK is Boltzmann’s constant. The first two terms in 
braces in equation (2.22) together comprise the familiar 
definition of entropy for fermions and bosons as reported by 
Nordheim [29]. The physical justification for these terms is 
that we must not only count particles but "holes" when we 
calculate the entropy [36]. The third term in braces in 
equation (2.22) is zero for fermions or bosons but for 
classical particles effectively subtracts the number flux 
from the usual definition of the entropy flux. Since the 


number flux is conserved, this term adds nothing to the 
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entropy production §& This term is also useful 


1 
Ala 


mathematically because when we choose 


2 ie ve 
ee 2 ney er 2h ae al ns ony 
we obtain EIN In (N,/A,) LOmPedlwtyoes. of pat iictes: 
Then, with Ne -k, equations (2.10) and (2.11) give us the 


entropy production equation: 


ae at — 4 ; 
S = -k 1 5 f [1n (N/A) IN,N,A,ARW d'V (BES) 
A Boltzmann H theorem can now be derived. We multiply 
equation (2.4) by NN, and integrate over both momentum 
variables. When we relabel some of our variables in the 


resulting expression, we obtain the following result: 


11 ¢s ik line ; 2.96 
5 | (SERN UM Ae ENE 0) (2520) 
Ly. SO kx 

where we have set = = (NAN, / AAR) : (NN, / AB) . We now add 


equation (2.26) to equation (2.25); consequently we obtain 
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We note that Neo Qmelsojs yo20) TS sequined by ithe exclusion 
principle for fermions, and 1s obvious for classical 
particles and bosons [12]. Therefore =*0 and hence 


In(Z) - £ + 150. Consequently, equation (2.27) tells us that 
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Ss") 4 70 which is Boltzmann's H theorem. When the entropy 
production is zero, Sa 7 we say that the gas is in 


equilibrium. 


C. General Decompositions 


Let us consider an arbitrary time-like vector field 


u(x") ; gut = 


derivative of u° by [31] 


= ss We can decompose the covariant 


Pe at genre eh } (2520) 
O a 3 


i & eee , 

where fae Cae GR Sei is the volume expansion, 
= jy ey AT ; is Nes ; 
Oe AS An - 34 y 04 ee) is the shear, Wa Hs Sellen is 
the vorticity, and a ype 0s the acceleration. 
In general, we can decompose the mass flux mM, with 

respect to u” by 

Oars: Qa BO a 

i = Gs ar ay (2329) 
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The quantity ni is called the mass density and the vector 


Ve is called the particle drift. 
Similarily, decompostion of the energy-momentum tensor 


with respect to u gives us that 
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The energy density is Pros the momentum flux is hh , the 
viscous stresses are rok , and we shall cal] P the bulk 


A A 


pressure. This quantity is the thermodynamic pressure Py 
Onl Vantn eau Mibr aun ot 1 dale 

The decomposition of the mass flux and the 
energy-momentum tensor of species A presented above are 
convenient for the analysis of the equilibrium and 


non-equilibrium states of the gas. We shall consider 


equi librium, farst.. it will besdisecussed gnychapter Li! . 
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In Kinetic theory, equilibrium maintained via 
collisions has a a special role; this role arises because 
the distribution function may be specified exactly. The 
specification of the equilibrium distribution function 
allows us to calculate the mass flux, the energy-momentum 
tensor, etc. exactly. These quantities and their balance 
equations describe the equilibrium behaviour of the gas and 
provide a foundation for the investigation of the physics of 
a gas in a state close to equilibrium. 

In this chapter we shall discuss equilibrium. We shall 
first specify the equilibrium distribution function and 
define the chemical potentials and the temperature. 
Secondly, we shall examine the restrictions on the gradients 
of the thermal potentials and the temperature imposed by the 
Boltzmann equation itself. We shall define some standard 
integrals of the equilibrium distribution function which 
allow us to calculate the mass flux, the energy-momentum 
tensor etc., immediately. This analysis provides us with the 
functional forms for the mass densities, energy densities, 
partial pressures, and entropy in terms of the thermal 
potentials and the temperature. Consequently we may obtain a 
Gibbs’ relation in equilibrium. Finally we obtain some 
expressions for the heat capacities at constant pressure and 


volume. 
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AD Thesequitior vm Dastribpution Function 

We have defined equilibrium to be that state for which 
the entropy production is zero. Let us denote the 
equilibrium situation by a superscript 0. When we examine 
equation (2.27) we note that zero entropy production 
requires that DEEDS On. Jhus In (N/A) is a collisional 
invariant. The most general form for In(N,/A,) is given by 


[33] 
In (N/A) = a,(x) + B(x) Pa, (3.1) 


Equation (3.1) now implies that 


° Ba ° g,exp(-c,-B P44) 
Ncsocomanns as abana: LS Ge aera ear ae f (3552) 


a = 
exp(-a,-8 p,,)-e, exp(-o,-B py, )-e, 


The coefficients a, and B, have physical 


A 
interpretations. The coefficient oN is the relativistic 
chemical potential and is related to the classical chemical 


potential u, by [14] 


a, = (m,/T) ( 1 + Ete e(Ce ease tun lis): (eo) 


To avoid confusion between the two types of chemical 


potential we shall refer to a, as the thermal potential of 


A 


species A. 
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We define a temperature T(x) by B= 1/kT where g=(-B 
We must have that 8° be time-like so that integrals of the 


distribution function converge. We also define a unit flow 


vector u~ by u = B°/B +; hence B = Bu’. The unit flow vector 


u” may also be considered to be the velocity of an observer 
moving with the fluid (comoving observer). For convenience 


we also define By = m,8 and Bn = mB". 
In equilibrium, the form of the distribution function 


implies that D 0. The Boltzmann equation now places 


N. = 
‘Some st 


Bes thice ions eon a, and Bs: 
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This must hold true for arbitrary time-like or null momenta. 


Hence =0 for all species. If the gas consists of only 


ALR 
massless particles we conclude that g° is conformal Killing: 


Re th an : 
= . However, if at least one component of our 
Boulvy ~ Dv] p 


gas consists of massive particles, then ge is a Killing 


vector: B =0. Thus the space-time is stationary. 


Cu] v) 
Since g is a Killing vector we conclude, after some 
a 


algebraic manipulation that B- Bu, Oe 0; anG* 4-0. hence 


the gas is also shear free and expansionless. 


B. Standard Integrals and Functions 
In equilibrium, all moments of the distribution 
function can be expanded in terms of a set of standard 


integrals of the equilibrium distribution function. These 
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standard integrals are defined by 


Qy.--%, — 1 Tp On, 
Th = el | NaPa Bie ks 
oA 
(355) 
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These integrals single out only one direction in space-time, 
u, + SO that they must be constructed out of scalars, the 
metric tensor, and the flow vector u- . Therefore, we define 


the following quantities [13]: 
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where 2q<n. We then have the following orthogonality 
relation: 
le Oey) Sanna tee eae (3.8) 
Meaty (q)a]---0, eerie coped 
Expanding equation (3.5) gives us 
ay On Bie G1. 0-G (3.9) 
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Contracting equa bions 3. cmamcedva. 0.) wititem io soueeciid 
employing the orthogonality relations (3.8) produces 
ae (3.44) 
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(22) 


Integral representations for these coefficients are deduced 
in Appendix A. Appendix A also includes general differential 
and recursion relations between these coefficients which are 


extremely useful. 


Ce the Physics of the Equa librium 
We can now examine the equilibrium form of the moments 


introduced in chapter II. In equilibrium, n, =N,, so that 


for the mass flux we obtain ai ges When we apply the 


decomposition (2.29) with respect to the flow vector u_ , we 


conclude that n, = Iyy9 and 5,20 ; furthermore, the number 


density becomes mn, = Dayo! Ma" 


Similarily, for the energy-momentum tensor we obtain 


HV Gem Ve S uv uv : : ; 
T, Ly Trot + T5434 . Comparing this result with the 
Qa 


decomposition (2.31) with respect to the flow vector u , we 
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the bulk stress to be 7,/3=P,-I,,, SO that in equilibrium 7,70. 
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Employing this expression in equation (2.22) gives us the 


entropy flux in equilibrium: 
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Israel and Stewart [16] have shown that Spa Wee lines The 
entropy is defined by Sy ==, 3, Noting that 
(e2 - 1)(n, ~ BP,)=0 for all species we then have that 

s,= ef oa, - Bee, +2, } (3.16) 
Defining O,= ka,/m, we deduce that 
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When we differentiate equation (3.16) and apply relationship 
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TdS, = do, - TO,dn, : ope hey 
The differentials appearing in this and subsequent equations 
in this chapter are constrained to virtual displacements 
between equilibrium states only. Multiplying equation (3.18) 
by thewarbitrary parameter veand Using» (3.17) Gives us the 


Gibbs. relation 
Td(S,V) = d(p,V) enh) TO,d(n,V) 7 (Bato) 


When we perform a Legendre transformation H,= (0,+P,)/V we 


obtain an alternative form of equation (3.19): 
HERG ID Gis AN adie) aaa) 3i2 OF) 


Let us now choose V=1/n where n= oe We then define 
the specific entropy per unit mass by o= oe , the enthalpy 
sive 1b = Ores , the internal energy by ie $P,/n eect 
ny, =n,/n whichis athe fractional propor lion Dy mass of 
species A. Employing these definitions and summing equations 
(3.19) and (3.20) over all species A produces the following 


forms tor the Gibbs relation: 
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For constant chemical composition ane 30m This requires that 
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where ne is defined by equation (A22). Summing this equation 
over all species and rearranging the result gives us the 


relation 
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FOP a Gas Consisting only o1 classical particles, 
i 230i ; sO we obtain, after some algebra, that 
ae = 0. Hence, constant pressure implies that the 
partial pressures are constant for a gas of classical 
particles. 


When we employ equation (A19) for I along with 


A20 
equation (3.22) we obtain 


J 1 a ¥ m,D 
(dp)- = - ) So Se “: EY ) ee | dT ; (ame 8) 
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where Da00 is defined by equation (A21). Inserting equation 


(3.24) into (3.28) gives us an alternative form of equation 
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The specific heat at constant volume is now defined and 


given by 
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The adiabatic index y is defined by 
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After some algebra we obtain the familiar result that 


V6a, (3.34) 

ine calculations. 1 fom equal 1omMs9 43.22) slo \o.o4) “Only 
apply to a system not containing photons. This is because, 
as the temperature increases, the number of photons in a 
fixed volume increases. Hence the proportion by mass of 
photons increases, that is, dn ,# OF FOR PROLOtSw nus sr Old 
gas of matter and photons, we apply the above calculations 
to the matter part of the gas and employ the procedure 
reported by Chandrasekhar [6] to compute specific heats. 

We now have a number of useful formulae which enable us 


to examine the non-equilibrium case, which is treated in 
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IV. Non-equilibrium 

The non-equilibrium state is much more complex than 
equilibrium and discussion of it will take the remainder of 
this thesis. In this chapter we wish to discuss the approach 
we will employ, namely, the fourteen moment relativistic 
Grad method. 

We shall begin our analysis by specifying the 
distribution function in the fourteen moment approximation. 
This distribution function will then be employed to find the 
mass flux and energy-momentum tensor in terms of the 
deviations from equilibrium which appear in the distribution 
function. However, the mass flux and energy-momentum tensor 
may be decomposed in terms of mass density, energy density, 
heat flux, etc. Hence we can specify the relationship 
between these physical quantities and the deviations of the 
distribution function from equilibrium. This relationship 1s 
then used to find the structures of the double-momentum flux 
and the entropy flux in terms of thermodynamic functions and 
the physical quantities appearing in the mass flux and 
energy-momentum tensor. Our non-equilibrium solution permits 
us to compare the actual distribution function with an 
arbitrarily chosen nearby equilibrium distribution function. 
Thus, in the final section of this chapter we are motivated 
toesinow that our descriplion’ Of =the gas iss independent ot 


this cnoice. 


4} 


77> 


‘ | eee: 
as soe 
7 
7 opitaahn in crusty Ee ca 
i ; vr 
4) sel omen ent wise Tae vate huts 
* a 
; i 1s ; A) “ 
ae) , art He ‘7 ij 7 j ul et 9 eyarh rie ait ne a5, 
7 ! = : _ 7 7 7 
Pa 
a A a , SS nf ort Fy ean vi: * +e J CATS air? ad ini a yaip 
_ ; sé - * i : 
’ Troe Aesscuet of? OVhemen: ..yol ane itw ow 
A 
revi tom 510 
wt onrv> los Va 2 sua suo Poet Uf grta ot 
; j 
oS # 5 @ th 
+ rT iv Y 4a Pp] Ay ww wold rie noe tad 9 nos 
i i. ti ie f 
: poe a% « 
yy bel? of beesl aes sai\tony late Mote oni ies yyetts ib @ bi 
« °¢ j j 
anes th Wen yvouuss bee AUD? zen 
j > a“ 2 = 
rik yr 4 4 - i ‘i. * er { a 1 i] i | ra ‘ ¢ ent Taine ab, 
4 . j= { ‘ee ‘. ri ij J | ’ .) Al cri fonu? 
Wiles br? ’ ¥ » es ve Goce wt 7 1 eit jes; iT my LULA al : 
j amy} gy \* att “WF f & = rl cist 7 
wt? Vo mi mn dak add ont biteeia Psarsyen VV : 
- 
; fey 2h met 7 mort ration? noiIuart af 
? ¢ Ye ‘ ' ap . 
. , . 
iT eee ai Ol ag14 ve j wis LNT. T. Ge beau Aer 
eee ale * . “a 
t si HUT OPO Ad ey) eee? nf eula vetting se brie 
: ° : an - 
oO ? 3eemn ny tii (esr) elt arias) Vey avi et! 
ei. 
Vo et nme j 3 ‘ arth #4 vi mOen 49 men ten ida i arr 
7 oe 
ne ht te tant AvP oid h3thn .eurae SI pene 
cet dort qui oda veth aus mdi | iupes vdben ek a" 


- 


balove tan ane SW relysno dirt to ootlimed I SATo oily 


bi , 
Jo rat a ae [ 
1, _ : 


r 


Pe 7. 


a 


, 


2n0 e)) 
7 


— ease 
Yo rors - paw ye a tsar we 


ay a a 


is 


. 


. _ 


42 


A. The Grad Fourteen Moment Method 

The equilibrium state does not manifest some 
interesting physical phenomena such as shear viscosity, bulk 
viscosity, and heat flux. To describe such phenomena we must 
discuss the non-equilibrium state. 

To examine non-equilibrium, we employ the relativistic 


Grad methoa of moments. We let 


y,,) = In(N,/A,) = In(W,/A,) + £, 5 (4.1) 


where InN, /A,) is given by equation (3.1) with g* (common 
to all species) replaced by B) (unique for each of the 
| species). We also have B. = But , uu, 3= ol and 8, = m,8,. The 
quantity eR 1s assumed to be small compared to In(N,/A,) and 
iLomodicCmecOu Dera iis .Oncelh (0,) 3 tat S tht Onms us tna we rare 
close} tO equi iipraum: 

The relativistic Grad fourteen moment approximation is 


obtained by assuming that a ey cl elUicleleene te pelalerg ela) (ong 


momentum: 


A ATP 
aN = a, Gs) + b, Go) wyy + cy (x)W a Mae CAS) 
atte ~X TAG ; 
The quantities a, » bys Cy are unknown first order 
variables which describe the deviation from equilibrium. We 
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To facilitate computation, we define the following 


Quantities: 


_ ae = 
b = ° = . 
ae en ee ge eee 
Eee ee Se ee (4.3) 
A DATA. 2 Ae Aah ATA ? 
= ON iL 
Shag Pen Mak 3 ee 5 
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- and c, are all orthogonal) to uP Sewer ch 


Qa 
Note that 2 » Cy 
these definitions we may rewrite equations (4.2) and (4.1) 
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(4.4) 
B 
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When we examine equation (4.5) we note that it is 
invariant under two classes of transformations of order one. 


The first class of transformations are frame changes: 
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Since the individual species are close to a common 


equilibrium, we expect that all of the u's are close to a 


A 
common u, that is u,= utd) where 6 is 0, . We therefore 
A A ig 
employ the frame changes, equation (4.6), to adjust 


equations (4.4) and (4.5) so that we have a common u° for 
all species. Once we have done this however, we still have 


the freedom of frame changes 


applied to all species. 

These first order transformations do not carry any 
physical content. The specification of an equilibrium 
dis tha butions UnecuLOon eG lose stontne actual sai strabucion 
function is not unique. The first order transformations 
allow Us to, transform toa different choice of the 
Cou word uM das (mt button. LUNnGhLOn. wuUneChOVece sofa thts 
function, which we "match" to the actual one, is therefore 
just a matter of mathematical convenience. 

We define the deviation from equilibrium, oN, , of the 
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Then, when we regard equation (4.1) as a Taylor series 


expansion of Im(N,/A,) about the equilibrium value 


In, /A,) in terms of oN, to first order, we have 
oO ° Cha 
In(N,/A,) = InW,/A,) Pe oN, : ee aie) 
. N.=N 
aN 


Hence we can immediately identify ON, : 
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B. The Non-equilibrium Tensor Structures 

We now have all the information necessary to deduce the 
expressions for the mass flux, energy-momentum tensor, etc., 
to. his UOrcder: 

From tnescderinitionm for mass. flux, equation (2.14 }" and 
employing equations (4.10), (4.12), the detailed structure 
Of 4,5 and the definitions (3.5), we obtain the following 


intermediate result: 
M, = 1, + a,J Ran ne Bae (A. 44) 


Now, using the definitions (4.3) and the detailed structure 
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Shecetatyons ino Mand (sat 0)er and comoaringe tie results 
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oe 2 + 6n : (45515) 


Qa a 
Urey Zunes (4516) 
where én, is given by 


b CA?) 


a= Jarota * JazoPa * Sago + Jaz oa 


We may split the mass flux into zeroth and first order 


Dames: 
MO = ME + MS (418) 
where 
Maa (4.19) 
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intermediate result: 


ete = el eaten eee os +e J , UO sat) 


Expendindeedlat tony 4221) wie cocualyonse(340) sand (eo. 10), 
and comparing the result with equation (2.31) allows us to 


obtain. the following results: 
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We split the energy-momentum tensor into zeroth and 


first order parts by letting 
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where 


°aB me a B ak i 
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Let us define the heat flux qa by 
is Ses er ee aa (4.30) 
This quantity is invariant to first order under frame 
changes, equation (4.9), (Israel and Stewart[16]). Then, 


from equations (4.16) and (4.25) we have 
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then we can find a, and By We call these the fitting 
conditions. This justifies calling the transformations 
(4.7), (4.8) fitting changes. We choose to let our fitting 
conditions remain arbitrary so that our subsequent analysis 


remains independent of any choice of the fitting conditions. 


C. The Solution for the Deviations from Equilibrium 
For massive particles, equations (4.17), (4.23), 
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dimensional linear algebraic system of equations which is 
conveniently expressible as a six dimensional 


transformation: 
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Therefore, we only have thirteen equations relating the 


a a ap : 
ae ye Pe Cyr oy ) to the thirteen 


Qa a we 


four teen variables (a,, b 
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. Hence the problem is not 
solvable with the given information. We shall have to employ 
a different approach to the massless case. Hence, we shall 
defer discussion of massless particles to chapter VI. It 
will now be understood that, in what follows, we are dealing 
solely with massive particles for all species. 

Let us denote the inverse matrix of the matrix 
appearing in equation (4.33) by (2) with matrix elements 
ee (i,j = 1°6). Then the inverse relation to equation (4.33) 
is given by 
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The structures of the mass flux and energy-momentum 
tensor have completely determined the variables appearing in 
Ee Therefore, the structures of all other moments of the 
distribution function are completely determined by the 
structures of the mass flux and the energy-momentum tensor. 


This is the hydrodynamical description of the gas as 


described by Israel and Stewart [16]. 


D. The Double Momentum and Entropy Fluxes 
We shall now discuss the two other moments which are of 
physical interest. These are the double momentum flux and 


the entropy flux, Let Us discuss the structure of athe double 


momentum flux first. 
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(4.49) 


As before, we use our Knowledge of the structure of Ny, Glin 
a to deduce the intermediate result 
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60. = én, + oc : (a5 2) 
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entropy flux to second order to properly discuss the 
transient thermodynamics of the gas (Israel and Stewart 
PG) (his willbe clarititedsim chapter Vv. 

Now consider the function 2,9 as given by equation 
(2.24) which we expand in a Taylor series to second order in 
ON ,: 


NG 


e ° °o o iL IN 2 
DR eel D1 arg 1s) NON is are oN) ae (4.65) 
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Evaluating this expression further gives us 
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= €,8,in{1-e,exp(a,+8,w,.)} as 2 f, : 


which reduces to the expression (3.13) in equilibrium. 


Inserting this expression into equation (2.22) gives us 
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Introducing the expressions for the mass flux, the 
energy-momentum tensor, and the double-momentum tensor we 


obtain an alternative form for the entropy flux: 
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kB, Cay if Pr) ~ On, : (4.70) 


The problem of deducing the structure of the entropy 
Flux now reduces to deducing the structure of On . We employ 
expression (4.2) to replace just one of the f,'s appearing 
in the right-hand side of equation (4.68). We obtain, after 


some algebra, that 
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To help re-express equation (4.72) we note the 


Following three relationships: 
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We now use equations (4.34), (4.55) and the three 


relationships immediately above to re-express equation 


(4.72) and thereby obtain Q» Tots inal Orme: 
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gr = S f+ , (4.88) 
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ao : zl 4.89 
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EE invariance: to First Urder 

We have deduced the non-equilibrium structures of the 
mass flux, the energy-momentum tensor, the double-momentum 
tensor, and the entropy flux. These four tensors describe 
alivofsthev physics of our mul tiscomponent gasm ihe 
description of the physics must, however, be independent of 
how we separate these tensors into zeroth and first order 
parts, that is, the description must be invariant under the 
first order frame and fitting changes. 

Let n represent the number of species in our gas. 
Suppose we are given dn sets of (measured) data (n,» P,) We 


select On numbers (a,, Bd and a rest frame u such that the 
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Tibtino cond) tions spectry the 2n numbers (én,, 6©,) and our 
0 


5 a 3 3 
choice of u specifies Gena 


The fitting changes and the frame changes tell us that 
our selection of (a,, B,) is not unique to first order. We 
could have specified another set (a,> Br) and another u'” 
since these are different from our first set by order one, 
lel US set 
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Then, invariance of In(N,/A,) requires that 
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Let us put equations (4.93) into Nt . We obtain 


Nae +AN where AN, = (NA, /8,) SE y+ Similarly, we can show 
that expression (4.12) becomes 6N, = ON,-AN,. Thus 
Ne ON a= N,tON, is invariant under fitting and frame 
changes. This decomposition of Ny is precisely the one used 
to obtain expressions for the mass flux, the energy-momentum 
tensor, and the double-momentum tensor. Hence, these tensors 
are invariant under fitting and frame changes. 

The case of the entropy flux is more subtle because it 


is evaluated to second order. Equation (4.65) may be written 


in the primed system as 
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Again, we expand the primed variables in terms of the 


unprimed variables via equation (4.93). Then we have that 
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Addition of the last three equations above now tell us that 


OR Og SON gO Mea rete at Uf aly Cro : (4.97) 
Therefore the entropy flux is invariant to second order 
under fitting and frame changes. 

We conclude that our description of the physics is 
independent of our choice of Rothe LOeiins t order: 
Consequently, we can choose any set of these quantities 
which is most convenient. However, WE note that it is 


impossible to choose, in general, a frame such that all of 
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's or the hats are zero. 
tHhereararne, two fitting choices that are most convenient. 


The first choice is 6n,= 6p,=0 for all species. This means 
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that Tai and Tyog are the actual number density , and 


energy density p, respectively. This has the advantage that 


A 
om is greatly simplified. However, we obtain a different 
temperature Gorveachs component — from this first 71 tting 
choice we may obtain the second choice by fitting changes to 
obtain a common B for all species, which means that we wil] 
have a common temperature for all species. Let (n> By) be 


the values which imply that én,= 89, =0 . Employing the first 


two equations in (4.93) tells us that 


iy lengthy hogy - 
(4.98) 
ae 0 eo 
We wish to maintain én,=0 . This implies that 
ba, = (Ty990/4Fa10) SB a so that Spy = - OD p99/Fa10) OB a . We specify 
B by requiring that 
rhe a areas Sid ) 6e, = 0 (4.99) 
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Since 68, = m,(B-B,) | equation (4.99) implies that 


2 { L@ 4208) Faro : { LO. 90%)/Tar0 } (4.100) 


This completes our discussion of the non-equilibrium 
structures and the invariance of the physics. We shall now 
turn our attention to the discussion of the transient 
thermodynamics of the gas for massive species. This is 


performed in chapter V. 
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V. Transient Thermodynamics 

sO far, in our discussion of the gas, we have not 
considered the detailed structure of the covariant 
derivatives of the mass flux, the energy-momentum tensor, 
the double-momentum flux, or the entropy flux. A full 
discussion of those structures will require us to compute 
covariant derivatives of the first order quantities 
Coe! Oe 


TT 


BUN as 


(On,, SP,» Ty, The retention of these 
derivatives in our calculations is a feature of transient 
thermodynamics as opposed to quasi-stationary theory. 
Transient theory versus quasi-stationary theory may be 
schematically explained in terms of length scales [15,34]. 
Let A represent the mean free path, the average distance a 
particle travels before it suffers a collision. Let L be the 
characteristic length over which macroscopic variables such 
as the number and energy densities change. We expect that, 
schematically, 1/Lvp'Y,N/N . Now let L’ be the characteristic 
distance over which first order quantities such as én, and 


69, change appreciably. We expect that 1/L' vp'V £/£ Pare 


A 
relative importance of these scales may be deduced from a 
schematic analysis of the Boltzmann equation. For a single 
Spcelles aswOteclassicdl particles, tnewecol ii sioOnetenm iiiathe 
Boltzmann equation is D NvN-of VNE/A where o is the 

could 
collision cross section. The Boltzmann equation is then 


schematically represented as piv Nw NE/A With the 


substitution N = N(I4£) for classical particles we have 
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phy N/N) Gs pi £/£)£ + phy N/N)E igi yf) eae Daa) 


In terms of our length scales this may be written as 

L/L + £/L“+ £/L > inis shows sihat fw X/i <<) thus first 
order quantities are about the size of the ratio of the mean 
free path and the characteristic length scale of macroscopic 
quantities. the third: term 1s of the size Vai which 71s much 
smaller than order one and negligible. Now if L’ jis about 
the same size as L then the second term is also negligible; 
this is the realm of quasi-stationary theory. On the other 
hand, if L’ is about the size of the mean free path, then 
the second term is also of size 1/L and is not negligible; 
this is the realm of transient theory. In summary therefore, 
when first order quantities change over a scale which is 
large compared to the mean free path then they are 
negligible and we are in the realm of quasi-stationary 
theory. However, if they change over a scale which is 
comparable to the mean free path then they are not 


negligible and we are in the realm of transient theory. 


A. Preliminary Discussion 
Let us begin our discussion by defining some 


terminology which will be useful later on. If we let woe 


represent any tensor function where a(n) = Os 20, then the 


time derivative and spatial derivative of ee are given by 
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Lneparteicu larg «id ye iS a Scalar efulactrom. of a, and Ba 
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Now on and By are defined in non-equilibrium by the choice 


Sof thercompoarison equilibrium distribution functions in 


equilibrium we noted that a = 0, B Ory; and 


A(u|v) ~ 
In a non-equilibrium state which is close to 


Alu 
eal BAU, j 
equilibrium we therefore expect that oN Veg 04> Nea pee = One 


and 8, = 8u+0,.. Hence equation (5.3) may be rewritten as 
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The time derivative of X, may now be seen to be of order 
one: lO B,)=0,- This property, along with definition 
(5.5) and equation (5.4) will simplify the computation of 


the derivative of the double-momentum tensor. 
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Felus also define the quantity iN by 


EeteuSswalsomnote a number oreuseful relationships: 


hy = huyu + AY bh, : haw, = - AAU, : (5iee/G) 
Oh (C40 

qa = G4,U + A neue Oe a (5,8) 
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We define the symmetric, spatial, trace free part of an 
ap 
arbitrary second rank tensor X,_ by 
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topics of this chapter. 


B. The Derivatives of the Physical Tensors 

InesTairstsand= zeroth order parts of athe: mass 1ux as 
given by equations (4.19) and (4.20) may be differentiated 
and contracted to give us the following two equations, 


Correct to order one: 
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We may solve for an in terms of hk 
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differentiated and contracted to produce expressions for 
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order: 
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In equations (5.24) and (5.25) we have retained 6c, and 
V,955 for the sake of simplicity. We can obtain expressions 
for these quantities by differentiating equation (4.55). We 
then have, correct to first order, the following two 


results: 
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Straightforward differentiation and contraction and using 


definition (5.2) provides us with the following result: 


oe . QM : Ce y Qh is ater 
5 Q,7,,00, 2s; qm, 60, 3 0,00 ,5P, oo 9,767, 
+5 Qin,ba, +> Qe6p,6n, + F Qede,on, + Q.én, On, 
a a Tr adn ii a TTI os TU, 
| a Ts fel yi oH TH is a Thaw, a 
rE ar Tce or aac i ae Taig 
s on vane a Ue ita e: ay Hy nth 
‘i a TP ala ‘3 a CN ‘ a" SPATa4y 


LES) wae UBS a -15 a 
=p HOPS N pear 0 URE een cra 


We can, however, compute Ant bY another method. If we 


differentiate the original definition of Ow equation 
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CV as ae 2 0. . 

ata | (dyeod} aC UN dv, ; (5.35) 


where we have brought the covariant derivative inside the 
integral. The integrand may be expressed as 


te} o fe) oO 


NA 
ne ee Cae AWA Oo A A ee 
CN Ae/Be) 1 1o%A i. 264 By AP { on J atata - (5.36) 


x ae i 
(9%. 00 advfmupe Jebfadee ile 
ay & aks, _ 
5 — 333 
Afzu 1 5 TO Se' WV 160 ns nAeregt ee x te ; . 
: _ bch 
‘ j 1 rw a5 Peat ts Pf rit j w a | rs 3) et 7 a) 


ao +. Say a i 
‘ rey va 


D } 
' i” ih ¢ ,™ wa 
' j , 
' 
! bar ar 
I nin eo Ce pe 
i A . . 4, R { 
a | 
Pt ‘ ihe ; af - { 
’ + \ 
| 7 “ A ne 
fh '\ . 
» > | : ; i 
i ee | F t. A F 
' ' 4 / j ie 
7 hy 2 ha f Cj 
; + : ‘ ‘ i , ’ 
f fj , 
ee | yn 
4 ae 
- ' i ‘ Aid a4 \ pe j 
‘ ‘ 
rv 
i C if + i « ‘ lA : f 
’ f 


‘ ’ } x { 8 ej (rE a | 
= : 
. To " 
- ee | ; er biee Fabo (larirerae ait + Tos tT it 


if ( ~ ios F 7 
ee CET PA Th Il ¥y ari niaice ow, (88,8) 


a ‘ | 
+. ~ _ 
a. } . ai ’ t i ' / { iy j \ a - 


sie" yt art it pwhas a8) drei woNGS erty, diet 7 
2s Bae anus acs ‘yer ely sat 


a 


73 


Tiesseconca term on the might mandeside of this ase cauation 
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(5.40) 


We now substitute into equation (5.40) equations (5.16), 
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the remaining variables via equation (4.34). The resulting 
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agree with the right hand side of equation (5.34). A 
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Perrorming the Gitterentiation one the tirs.t. term ian. this 


equation eventually leads to the following result: 
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To check this expression, we compute the entropy production 
by another method. From the definition of the entropy flux 


we have 
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We may now substitute into equation (5.55) equations (5.13), 
Porm re Ont) be O01) same ee See Omari (Seo eae Well! 
as employing equation (4.34) to replace the remaining 
variables. This whole procedure produces equation (5.51). 
Hence we confirm that equation (5.51) is the correct 
expression for the entropy production. 

Equation (5.51) may be re-expressed for better clarity. 
We add equations (5.18) and (5.20) and rearrange the result 


to obtain the following expression: 
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In the case of the single species gas, Th 1a 0, 
equation (5.57) reduces to the result as given by Israel and 
Stewart [16]. Also, when we neglect the transient terms, we 
obtain precisely the stationary result in kinetic theory 


given by Israel [12]. 


D. Transport Equation Preliminaries 
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This 1s one relation that we shall require. The other 
relation that we require is obtained by noting that the 


thermodynamical relationship (Israel [12]) 
(aa) 


may be rewritten by expressing the differentials as spatial 


derivatives: 
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We now substitute equations (5.70) and (5.72) directly into 
the right-hand side of equation (5.22). We eventually obtain 


the following expression: 
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We now have all the information we need, in the appropriate 


form, to compute the transport equations. 
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equations (5.69) and (5.24). We then have that 
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The coefficients introduced here (and later on in our 


so he indices 


discussion) are figuratively of the form ta 


1 and j mean that nee is the jth coefficient in the ith 
transport equation. 

The left hand side of (5.74) must be evaluated via 
equation (2.21). Since this procedure is rather complex, the 
evaluation of the right hand side of equation (2.21) is 
deferred to chapter VII. From chapter VII, therefore, the 


left-hand side of equation (5.74) may be rewritten via 


equations (7.58) and (7.56). We therefore have that 
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Let us now replace the left-hand side of equation (5.74) by 
the right-hand side of equation (5.82). This gives us our 


first transport equation: 
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F. The Second Transport Equation 
To obtain the second transport equation we add 
equations (5.73) and (5.25); we therefore have that 
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The left-hand side of equation (5.86) may be rewritten via 
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Then, replacing the left-hand side of equation (5.86) by the 
right-hand side of equation (5.96), we deduce the second 


transport equation: 
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The third transport equation is obtained by adding 
equations (5.23) and (5.26) which gives us an expression for 
Ce ay A Ee ; however, this quantity is also given 
by equation (7.60). Equating these two expressions presents 
us with the third transport equation: 
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We may formally solve for the bulk stress, the heat 
flux, and the viscous stresses by noting that the transport 
equations are matrix equations in the species indices A and 
B. Then, if we assume that xX,5; (S335 ~6) has a matrix 
inverse ie, COD , we multiply the transport equations 
by these inverse matrices to obtain the following three 


versions of the transport equations: 
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= al . 

hap = XapsJpovip/p 3 (5.104) 
et 

MAB ARG SRST (5.105) 


The three transport equations above reduce, for a 
single species gas, to the forms reported by Israel and 
Stewart [16] provided we apply the fitting conditions 
én, = 6p, =0 and choose u by setting hi=0. We also note that 
the transport equations are linear, first order, coupled 
differential equations which in the single species case are 
hyperbolic and lead to propagation speeds of viscous and 
thermal effects less than that of light (Stewart[34]). This 
suggests that for the multi-component gas the transport 
equations are hyperbolic and describe causal propagations of 
the heat fluxes and viscous stresses; however, a proof of 
this is beyond the scope of this thesis. Our discussion of 
the massive components of the gas is complete. We now turn 


our attention to the light-like components which we wil] 


discussmin chapters Vik 
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VI. Massless Particles 
Massless particles, such as photons or neutrinos, 


possess zero length four velocities, that isw. w =0 (we 


o) 
ee ab 
reserve the symbol L for massless particles). This condition 
immediately implies that the energy-momentum tensor for 
these particles, as defined by equation (2.16), is trace 
free. We have already noted in chapter IV that this result 
deprives us of one of the conditions necessary to obtain a 
unique solution of the non-equilibrium situation, for 
massless particles, in the Grad fourteen moment 
approximation; this is because the number of linear 
equations relating fourteen variables is reduced from 
fourteen to thirteen in the case of massless particles. 
Hence, we must modify the fourteen moment relativistic Grad 
method in an appropriate manner to obtain such a solution. 
In this chapter we shall discuss precisely how to 
modify the fourteen moment method to obtain the 
non-equilibrium solution for massless particles. To achieve 
this we shall present our discussion as a miniature version 
of the previous chapters on massive particles. We will start 
out by specifying how the Boltzmann equation must be 
modified to account for the interactions between massless 
and massive particles. We introduce a change of variables to 
facilitate subsequent calculations by expressing the 
massless particle’s four-momentum in terms of a frequency 
and a unit spatial vector, both of which are defined with 


respect to the flow vector u-. This change of variable 
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allows us to describe the physics in a frequency dependent 
format. The equilibrium distribution function for massless 
particles allows us to define frequency dependent integrals 
and thermodynamic functions which will be employed to 
analyse the non-equilibrium form of the distribution 
TUnCiLIGn.s Ihis distribution function, 15 calculatedrVvia a 
frequency dependent version of the relativistic Grad method; 
this is obtained by assuming that the coefficients in the 
power series expansion in the particle’s four-momentum 
depends on frequency as well as on position. The 
non-equilibrium form of the distribution function, when 
inserted into the Boltzmann equation, leads directly to the 
transport equations. Finally, we shall derive the entropy 
production and define the coefficients of bulk viscosity, 


shear viscosity, and heat flux. 


A. Creation and Annihilation Processes 

Let us begin our analysis by considering the Boltzmann 
equation for massless particles. The collision term, 
Dae ake stated in chapter II, only takes into account 
binary collisions. binary collisions, bysthemselves;edo mot 
adequately describe massless particles. For example, the 
number of photons in an equilibrium state is a function of 
the temperature, that is, the number of photons from state 
to state is not conserved. We therefore introduce into our 


description of the gas creation and annihilation processes 


in which the number of massless particles is not conserved 
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[1]. These processes are figuratively of the form 
AtB+L > A*+B* (annihilation or absorption) and A+B > A*+B¥*+L 
(creation). 

We describe annihilation and creation processes for 
massless particles by the transition probabilities 
WP psPy oP, [Pa ePp) and WP, sPp Pg sPp oP) respectively. Then, 
following our analysis of chapter II, the contribution to 


Doo, by these processes is given by 
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where div = av, dv,dV,dv,, Then the Boltzmann equation for 


massless particles is given by 
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where B(N, ) , the contribution due to binary collisions, is 
given by the right hand side of equation (2.5) for A=L and 
the sum is taken over all massive particles. We are assuming 
here that interactions between massless particles are 
negligible (e.g. photon-photon and photon-neutrino 
processes). 

For the sake of consistency, we should also add the 
contribution of creation and absorption processes to 
D ola for massive particles. This contribution is given 
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Then the Boltzmann equation for massive particles is 


N = B(N,) se TOV) ; (6.4) 


Aeon ceo 4 
where B(N,) “is given by equation (2.5). 

The master balance equation (2.7) retains its validity 
here. We just have to ensure that Daa inc ludes T(N,) or 
T(N,) . AS a consequence of these extra terms, the master 
balance equation implies that the massless particle number 
(mass) flux is not conserved but the number and mass fluxes 
for massive particles are conserved. Furthermore, the 
energy-momentum tensor of the whole gas is still conserved 
and the Boltzmann H theorem is still valid [12]. Then the 
equilibrium situation as described in chapter I] remains 
valid. The equilibrium restrictions on the thermal 
potentials on require a, =0 if we include creation and 
annhilation processes. Hence, for massless particles we 
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Spectral (frequency dependent) forms of the number flux 


and energy-momentum tensor may be defined as 


Ne (vy,) = | NW, ¥, do, : (6410) 
ap 2 a p 
Te ie) = | Nw, Ww, v7 dQ, . (6. 11} 


When we integrate over frequency vy, we just obtain the usual 


number flux and energy-momentum tensor: 
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the spectral case: 
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Then, integration over the frequency just produces equation 
(3.5). Applying the analysis of chapter III to this case 
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C. The Modified Grad Method and Non-equilibrium 

We may now consider the description of the deviations 
from equilibrium. We shall assume that the flow vector u- 
has been chosen via frame changes to be common to all 
species. We expand In(N,/A,) , as in the massive particle 


case, by 
In(N, /A, ) = In(N, /A, ) oe bos 13) 


where we assume that f is given by 
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assumption here is that the coefficients ars bb and — 


which appear here are also functions of frequency as well as 
of position, whereas in the massive case (Grad fourteen 
moment approximation) they were assumed to be functions of 
pesition alone. 
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expression of 11,2]. The number of unknown variables 
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Therefore, the deviation from equilibrium is completely 
specified by the spectral energy-momentum tensor. 

The spectral number flux may be evaluated in terms of 
equations (6.22) and (6.24); comparing the result with 


expression (6.13) gives us 
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which in turn requires, from equation (6.25), that 


én (vy, ) = on, (vy) - Jy 99 6%z,) 98, : owo) 
Integration of these two relationships over frequency and 
emo loving derinitions (6.26) produces the wsual ti tting 


relationships (4.98). Although we may choose BE such that 
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7 it~ 9, no choice of 8 can make ¢n,(v,) or 6p, (v,) 


equal to zero. 
For each massless component of the gas we define the 


radiation temperature T B. = 1/kT, to be given by the 


LR eR 

Fitting conditions én, =5p, =0. We consider the matter part of 
the gas to consist only of massive particles. Therefore we 
define the matter temperature a, BR, = kr, to be given by 
the Trilling conditions én ,=0, 790 ,=0, where the sum is taken 
over all massive particle species. This case was previously 
discussed in chapter IV. We also define the Eckart 
temperature oe B= 1/kT,, by the fitting conditions 
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These conditions are the generalization of the Eckart 
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D. The Spectral Transport Equations 

Let us now derive the transport equations for massless 
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conditions imposed by the Boltzmann equation on the 
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preliminaries. 


Let us differentiate oe and a as given by equation 


(6.5). This provides us with two very useful results 
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we and i fixed. If Xr, 1Sea_ TUNGLION Only oT position and 


frequency, then the last term in equation (6.38) is zero. In 
this case, diffentiation holding frequency and ky fixed has 


the useful property that 
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which produce upon integration over frequency, via 
relationship (6.39), the usual time and spatial derivative 
as defined by equation (5.2). 
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To obtain an expression for the left hand side of the 
Boltzmann equation (6.2) we differentiate (6.43) and 
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to compute the first two terms, it 1s helpful to Know 


that 


-3 a 
(vy en . 


1 wae 6 (3 
= { a) oe He (Gina) see Te Die » } (6.45) 


Vv 


With the help of this expression and equations (6.38), 
(Bet0erand) (6.44 )eetthe stipstebermalny ecuationm 6244) 1s 
given by 
(OL yg ae ae I. Gs SRG ene tee ple 
cae EY Ae, Abhi Sie ee ee 
ie 


+ Ke LV , (v) + [3p, (,) = v, Pp (v,) Ju} (6.46) 


NlR 


a, B 
kk (3p, ,) = Vp, cvs Baas 


eit to able ibned Mel set onume peg, T | 
baa (88,3) sia vanesstle So) io ae cnet 
oe atte th pei 


ia it penncyongot = Poplat aie 
ere #2 rigs r aa Sa at eg COR of 7 


(Bb. 3) 


| | 0 | | ao 
{sphere * pay! jretggth = + is 7 

¢ 7 ao “Ty i. 

1 LS sega * hf Ga whe e" ae 


Yd mie) (88,0) eearqxenet I isda ow jentww? beese1qg OTe aa 


ent ontt of enmotese gs nid iu2ss adi efdmezes nent bons mist 7 
(RD,3) mottsups tot Not easrqKe berteeb a 

wood of futqisn at tf ,emnot owt Serr? ett etugmoo oT rs 
tend aU 

at 7 

(eh .d) ; { (calor ae + (be) St +6 } a = aio 2? ert) : i 
a ee | 


(86.8) amotisups bas notzesigxe 2atAt to qferd ey As 
at (03,3) notisuips ni mast Javit edt .(t8.3) bas ) 
} (; i 


The second term becomes, with the aid of equation (6.45), 


1 
TOL yg ve Shy Oe iy Lo Deke 


iis 


MBs, 


Br . a, By, 45 : 
2 a pg Loh, (ug) + Kt yap De) f- 


i 


For the remaining terms in equation (6.44) we obtain 


r a Nes 
ae Nee Oe { i th)! 


(6. 
a, B , ¥ 5 
Ft OD ee why ae) ; 
: Sa ee i peen a| 
Orleans gal Ey ar (6 
r oO a, B x , 
‘ ci ws Pre oka - Soil ada} : 
Bev.) = yf Ket (J 
Seay i lee ‘ le teria. wn i 
Bay : 
a Tl Tag” L aay, = iy Tap %y)} a 
ap e 
ae eS _ “4 al TT Oy, yu, ] y 


g A Oanbsey 
~ gk! Fe, Ope esas ts sn Tap %y)} 


104 


47) 


48) 


.49) 


DU) 


Sh) 


We may now substitute equations (6.46) to (6.51) into 
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calculation to chapter VII, which is devoted to such 


calculations. We extract the result of that calculation from 


atl -; ips ati wis it nee weenie > 


aved ew 94 
= 


i yak 

» .. ‘ f cy 

ij I J M5 as 
' 


F Lice ‘ y— | we) Hf be | 
J rd i 
Lis vl 
. ;@ i 4 a = ; a int py ‘ab = 
* | WP t DA ; , 
ij 1 ) At 
! 
ca al 
‘ ! ! f 
(@! 
. h 
: ] —- is Py 
if y bf Pa rv { 


she ts 5 


he : 7 " 
eared atrgit Sts assed biol ii 
Te | ie a ’ 


106 


equation (7.109) and equate it with equation (6.52). We 
obtain 
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equilibrium, we multiply equation (6.57) by 1, see Ky 

Kr KY, and integrate over d2. . We deduce that the four 


L 
restrictions obtained by this process are given by 


De ah =e 2 (6558) 


A pe a Bee ee (6.59) 


Ds cag> = Ds cog> ; (6.60) 
(au, VA) th (aB,yU,VA) a { 
Deas A is Ds apye IN OOIN an ye enol) 


When we substitute in the appropriate expressions, from 
Equationss 6.53) Om tO 5Olm ancenhommetsiaO)e LOU sre 2) yamebic 


first three of these results become 


oy) + F130.) - vppfQ,)16 + mAO) .-v hy, 0" i 


: {M6 Seer) Fm an aan 8PM 3M} ; 


a shih 


mot, ark? LAIVSe etth 16 ants 
ei yn oh (va 4 mel 

_ att tard woukeat oi pee 

vai covig on nsec fem) NAH 


ial al | Acc ae 
| ; 
(on.a* «MG Le aha 


3.3) ; u 7 vaca 


| 7 Fé ayy ta i § (hy eds 
ir2.% i A's ue © h ‘7 


ayaryt onahene Wino e1s (Qaryae at itt abul i Tadie gw 
st). (2 ReGt tel el nds® tne RE ot (ba Banat 
ensoat tlacatd salad ths 


P= 9 


Ny) : 


Bo 


lee ‘ (vill 4° BEL, ag p+ (gO) 
evans ee Zz ey We 


asians rsa ate 


Oy 


410, (,)+ 4 S13, (v,)- -v, 07 (v,) ut om) (6.63) 
* a n . Lok . 
oe ale ON) ae "ho ae Zz etELO (vp du, 


ee a LOH aPag ene es 


4 
(36, O97) Ie Yea] p> 7 ca 587 Shea YL) p> 
eee ap cree (6.64) 
Seay ae oe ee Oe es “40 gs 


KT 3" Lape -2 BLN? 


The fourth result may be rearranged to produce 


ts ia % 
2 { Tae ae “Tapia a° Les OR 
2 S AV; 
ee {50 "ha Og - “Marte ve ee an uaa (6.65) 
, : Y a 
a2 {> age Ct ~ “Mavta} ae Wea ie 
We may recover the stationary case, with Wye O Thon 
equations (6.62) to (6.64) by neglecting the transient 
terms. Noting that 
3199) — YET20 LD) = B3t30O%D) > (6.66) 
we obtain for the stationary limit 
[J Benn Cae Dsl ie + [J} 49%, /3kT, 10 
(6.67) 
Se aes © L EN aro ey ect) at 
us 
a IS ae WV, fons ee A) 
(6.68) 


Ser to ue ae LO nana ee ahee ae 


By, 5 enaitit oy : vont te 
(28.4) Jae via as ih lg pee 
sage <iger be i ga yt ys = mains! ehss + a 


gervi, , o a Atty sens Sytena? rela ae Sale” Node bal 
jnadanert: ont ghktootour co '#8Oi of (83-6) anole ? 
seit! on rei em 


ra. 8! 0 «Reta stat = Ae pag acs a 
i : . 

ae 

frit yenopteta att tw nbeiees 


on BE FAL et gh neha a 
eae ant? Beat cae . 
Aire . | 
_ lg ao 


one i m sls a Ai + 4 
fe © ts Gee ae . : - 


108 


i 
31 Ip 39%) / A5kT, Ia g (6.69) 
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If we choose T=T, and neglect the deviations from 
equilibrium of the massive particles, then we obtain, for 


the case of photons, the results reported by Straumann [35]. 


E. Frequency Averaging 
To integrate the transport equations over frequency we 


must define some way to average k and Ma; over frequency. 
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In the stationary case for photons where the contributions 
by massive particles is negligible, we would just divide 
equations (6.67) to (6.69) by the appropriate Ky and 


integrate. The appropriate average in that case is the 


Rosseland mean [3,35]: 
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We shall carry this averaging technique over to the 
transient case. 
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Scalar cOetiicient. Forstnese reasons we Shall takeedl lor 
Our “absorptlom coefficients: to be Scalars. 


Let us define the following frequency averages: 
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These are our transport equations. In the stationary limit, 
with ee and neglecting the deviations from equilibrium 
of the massive species, the first three of these equations 
are identical with the results of Thomas [37] and Anderson 
and Spiegel [2] recast into our notation, provided we assume 


that all of the absorption coefficients are the same. We 


also note that the first two of these equations are 
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expressions for -u T and kT 
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F. Coefficients for Dissipative Processes 
To identify the coefficients of bulk viscosity, shear 
viscosity, and heat flux, we must examine the entropy flux 


and entropy production. The entropy flux as defined by 


equation (2.22) becomes 
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Since decomposition of oN) up to second order, equation 
(4.95), is invariant under fitting and frame changes, 
equations (4.96) and (4.97), then equation (6.78) is also 
invariant under these transformations. In this case, we have 
just reproduced the massive particle result for the entropy 
flux, equation (4.67), except that in this case O, 1s given 
by equation (6.79); however, this latter expression can not 
be credibly integrated over frequency to produce an 
expression like (4.77). 

Let us compute the entropy production. We simply 


differentiate and contract (7.78). This produces 
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When we substitute these results into (6.80) we obtain 
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To determine the bulk viscosity we must examine the 
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simplicity that the deviations from equilibrium for massive 
particles are negligible. We shall also take T, sta Drauss 
case, equation (5.54) for species A becomes 
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When we sum this result over all massive species we obtain 


the entropy production for the matter: 


Oe re ne 


Qa 
le = kB Ty, Er (6.88) 
Since the total energy-momentum tensor of the gas is 
conserved we must have that 
CU haamees, an 6.69 
Ty Ir ae A ( 


Hence, summing equation (6.83) over all massless species and 
adding it to equation (6.88) gives us the total entropy 


production of the gas: 
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since ae for a common temperature. When we substitute 


this expression into equation (6.90) and define 
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This expression may be more conveniently expressed. 
When we hold the individual number densities constant, i.e. 


dn, = 0 , we obtain for any thermodynamic change that 


da, = My (Tn50/In1 0) dB, (6297) 
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Let us now sum these expressions over all species in the 


gas. This provides us with two convenient results: 
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Weomay OW SUDSTitUte these two restilts) ifite expression 


(6.96); therefore, we obtain the following result: 
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We note that in the stationary limit and for only one 
massless component in the gas, we have precisely the result 


for the coefficient of bulk viscosity first derived by 
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Weinberg [41]. Hence our result is the generalization to the 
case where more than one massless species is in the gas. We 
note that if the matter is itself highly relativistic then 
the matter pressure is about one-third the matter density, 
16. Pe p/3; In this case the bulk viscosity vanishes: 

We noted previously that, for the sake of consistency, 
we must include the contribution of the interactions between 


massive and massless particles in p Consequently, 


see ee 
such contributions must show up in the transport equations 
for massive particles. This contribution is computed in 
chapter VII and leads to equations (7.78) and (7.79). These 
contributions have exactly the same structure as the 
contributions due to massive particles, equations (7.41) and 
(7.42). Thus, all we have to do is incorporate equations 
(7.70) and (7.71) into (7.41) and (7.42) by redefining the 
Lapi and Xap; by 
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Then the transport equations for massive particles retain 
their structure as reported in chapter V except that now 
they include interactions with massless particles. 
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Structures that we deferred in this chapter and in chapter 


V; these computations are performed in chapter VII. 


VII. Non-equilibrium Conservation Laws 


In chapters V and VI we required a knowledge of the 


apr 
A |r 


probabilities and the deviations from equilibrium. These 


nN ; ees 
structures oa and U Tne LeErmMs OFC Lhe Coll i1si0M 


A 
structures were employed to write out the final forms of the 
transport equations. In this chapter it will be our task to 
deduce what those structures are. 

Our calculations in this chapter are of a highly 
technical nature; thus we shall present a brief overview of 
how our calculations done. 

To perform our calculations we first specify some 
assumptions and perform some preliminary calculations to 
facilitate subsequent calculations. Then we shall split up 
our discussion into three stages. The first stage consists 
of examining the collision structures of 1° es 


A AI 
massive particles due to interactions between massive 


i and wv for 
particles alone. The second stage consists of calculating 
the additional terms necessary to include the interactions 
between massive and massless particles. The third stage 
consists of determining the structure of De for 
massless particles due to interactions between massless 
particles and massive particles, considering interactions 
with other massless particles to be negligible. 

In each of these stages cited above we will proceed by 
expressing the distribution function in the collision 
integrals via the appropriate relativistic Grad method. The 


resulting expressions in terms of the equilibrium 
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distribution function and the deviation from equilibrium, 

fy: permits us to define several classes of integrals to 
simplify the expressions. We analyse these integrals in 
terms of their irreducible structures. Finally, these 
irreducible structures are introduced into the collision 
integrals to obtain the final expressions required to obtain 


the transport equations in the previous chapters. 


A. Assumptions and Preliminary Calculations 
Let us begin by assuming that we have detailed 
balancing for all binary collisions, creation processes, and 


absorption processes: 
ie ae ae fe ) 


* * * 
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This assumption is justified by the microscopic 
reversibility of these processes (Anderson [1], Israel 
[12]). For convenience we shall denote the transition 
probablities by Wap and Wap In equilibrium, the 
assumption of detailed balancing and the functional form of 


implies that D Na0 SOR atl 
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species. Consequently, a, and UR , are identically zero 
in equilibrium. Therefore in non-equilibrium we should 
expect these quantities to be of order one; furthermore, we 
expect them to be expressible in terms of the deviations 
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Before we enter into the main discussion, let us derive 
some preliminary results. When we employ the most general 
fitting conditions, each species has its own thermal 


potential Oy and inverse temperature B, However, each 


is close to a common 86 to first order, that is, we may write 


B, = BtoB, where 6B, 18 order one. Therefore we conclude, 


after some algebra, that for binary collisions 
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where we have employed the relation B Tp, 1=0 which implies 


that the quantity which appears in the exponent in equation 


(7.2) is order one. Since [a]=0 equation (7.2) implies that 
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Let us now subtract equation (7.6) from (7.7), (7.8) 
from (7.9), and (7.10) from (7.11). We obtain for binary 
collisions, absorptions and creations respectively 
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We may express the Boltzmann collision terms for 
massive and massless species via equations (7.12), (7.13), 


and (7.14). Let us reserve the subscripts A, B for timelike 
particles and the subscript L to refer to massless particles 


only. Thus we have for massive particles 
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interactions with massless particles, is given by 
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Since our discussion now becomes more complex, we shal] 
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B. Case One: Massive Particles Alone 
We now consider equation (7.15) in which we now assume 


lakene 36 is negligible. This means that we are only 


ror 
consieaering elastic binary collisions of massive particles 


with other massive particles. If we multiply (7.15) by a 
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If we rewrite the right hand side of (7.18) by exchanging 
starred and unstarred variables, employ equation (7.2), and 


then add the result to (7 718)8we obtain 
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Equations: (7.19)) and! (7 20jenow suogest that, tor 
convenience, we should define an integral KE by 


a il! °o fe} CHOK + * 
ABCa(n) 6a) Se th aes Cay “ca 
- area e 


= oe J Ky 
W pe Ww hey FE Wap x 


This integral is symmetric in the indices a(n) and has the 
following properties: 
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These tensors have the following properties: 
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To proceed further we must deduce the irreducible 
structures of the four tensors which are defined by 
equations (7.28) to (7.31). These tensors must be 
constructed out of scalars, the metric tensor, and the flow 
vector u” . We note that our selection of a common flow 
vector for all species has greatly simplified the complexity 
of the tensor structures. We deduce the tensor structures by 
writing down the most general form that the tensors can have 
and then simplify by using the properties (7.32). 
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protons, and neutrons. Let us therefore assume that these 
coefficients are, like the number and energy densities, 
Known quantities. 
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We note that if we choose the fitting conditions to give us 


a common inverse temperature 8 then L via the 
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C. Case Two: Massive and Massless Particle Interactions 

The five equations immediately above do not include the 
interactions of the massive species A with massless 
particles. We now wish to Know how they must be modified if 
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As before, we are primarily interested in the two cases 
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D. Case Three: Massless Particles Alone 

The last case we have to consider is the evaluation of 
D 61. Let us reproduce equation (7.17) in a more 
convenient arrangement: 
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then equation (7.90) simpiarres to 
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To make further progress we must express Ko, and K 6), 


in terms of their irreducible structures. Since they are 


tensor functions of position, momentum and the flow vector 


Q ; 
u , we obtain 
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NLA2 ~~ thn a EE) } (7.96) 
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then we have for equation (7.94) the following result: 
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This quantity is zero if all of the inverse temperatures, 
Bos ar the same. 
We must now determine how to evaluate the last integral 


in equation (7.102). Since we are integrating over the 


variables Pa Pe Pr Wrand an ion trotetiesstructuresct 
f preguationy (6525) —. thistimntedralamus tbe wa erunctron 


ie 
solely of the variables Py and position. We conclude that 
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where the coefficients K;} must depend solely upon position 


and frequency. 
Equation (7.102) becomes in terms of (7.103) and 
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(6.24). Hence, when we define 
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solve for the massless particle case transport equations. 
In summary, we have developed in this chapter three 


sets of equations. The first set describes i and U of 


apr 
A ie 


the massive particle components excluding interactions with 
massless particles. The second set of equations describes 


the additional terms which we have to include in ite 


A 
aBA 
Ke 


is and 


U if we wish to include the reactions of massive 


particles with massless particles. Finally, the third set of 


equations describes p aces massless particles reacting 
co 


with massive particles, reactions with other massless 


particles being considered as negligible. We have also shown 


aA 


BA 
that the plausible expansion of ripe th Jas and D n 


aoe 
terms of the deviations from equilibrium is quite valid; 
furthermore, the calculation indicates how the scalar 


coefficients are related to the transition probabilities. 
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VIII. Summary 

In this thesis we have developed a non-stationary 
theory of mixtures within the context of Boltzmann kinetic 
theory. The two major results that we derived are, for each 
component of the gas, the expression for the entropy 
production, and the transport equations which describe the 
behaviour of the thermal and viscous effects. These results 
are important because the entropy production and the 
transport equations, in conjunction with the balance 
equations for the mass flux and the energy-momentum tensor, 
are just the equations we need to describe the macroscopic 
behaviour of the gas. 

The transport equations indicate the synergistic 
effects in the gas. Not only do the thermal fluxes and 
viscous stresses of a particular species depend upon purely 
thermodynamical quantities defined with respect to that 
species, but they also depend upon the thermodynamical 
quantities, heat fluxes, and viscous stresses of all the 
other species in the gas. In addition, the transport 
equations include relaxation terms to describe the transient 
effects in the gas. These relaxation terms also appear in 
the entropy flux and imply that, for any initial deviation 
from equilibrium, the entropy is reduced, 1.e. equilibrium 
is a local maximum. This is a manifestation of positive 
relaxation times. 

Another feature of the transport equations is their 


insensitivity to the nature of the collisional terms 
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appearing in them. The transport equations may be written 
schematically as TD=COLL where the lefthand side, TD 
(thermodynamical), contains thermodynamical functions, their 
derivatives, heat fluxes, and viscous stresses etc. The 
Prohthand side, COLL” (collisional). are the coldisiona! 
structures formed from a linear combination of the 
macroscopic deviations from equilibrium and coefficients 
determined from the collision cross-sections. While these 
collisional coefficients are sensitive to the nature of the 
collisions involved in the gas, the general algebraic 
StUCLUREeG O17 COLE Ws not. Furthermore.) cur mains result, tor 
both the massive and massless cases, is the specification of 
TD for all the transport equations. These structures are 
independent of the details appearing in COLL and hence are 
applicable to all situations whereas .COLL will change 
according to the system being considered. 

Our theory of a multi-component gas with transient 
effects is a generalization of, and an improvement over, the 
single species non-stationary theory and the 
quasi-stationary multi-component theory. Under the special 
COndiltvOns Of a Gas Consisting of only one massive particle 
species, our theory is identical to the non-stationary 
theory of Israel and Stewart [16]. On the other hand, when 
we neglect the transient terms in our theory we 
produce,assuming the vorticity is zero, w= 0, the stationary 
multi-component theory of Stewart [33]. Finally, for the 


massless case, when the deviations from equilibrium for the 
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matter may be neglected altogether, we obtain the 
quasi-stationary results of Straumann [35] when we neglect 
the transient terms in our theory. 

Our theory is also an improvement over the 
quasi-stationary multi-component case and the single 
component transient case because it has a wider range of 
applicability. The situations envisioned here are cases 
where a multicomponent approach is necessary and where a 
quasi-stationary theory is inappropriate. Possible 
applications for our theory which come to mind are of an 
astrophysical nature: the accretion of matter through the 
boundary between a star and a hypothetical neutron star or 
black hole imbedded in the star’s core (Thorne and Zytkow 
[38]); accretion disks around black holes and neutron stars; 
and the leptonic era in the early history of the universe. 

In’conclusion, therefore, we have a non-stationary 
theory of a multi-component gas containing massive and 
massless particles which is a significant improvement over 
previous theories; and this theory is applicable to 
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Appendix A: The Standard Thermodynamic Functions 
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Appendix B: List of Symbols 


-a deviation from equilibrium appearing in £, \p- 
oe 

-a numerical coefficient (p. 34 ). 

-a deviation from equilibrium appearing in ey er 
a3) 


-a deviation from equilibrium appearing in ¢£. (p. 


A 
San 

-a deviation from equilibrium appearing in ¢£, (p. 
GO. 


-the collisional terms due to binary collisions 


alone (pe 91 ye 


-the speed of light (=1 in this thesis })(p. 17 ). 


-a deviation from equilibrium appearing in ¢£, (p. 
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-a deviation from equilibrium appearing in ay 2: 
A320), 


-a deviation from equilibrium appearing in £. (p. 
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Ae le 

-a deviation From equilibrium appearing in £, (Pp. 
a2 

-the specific heat at constant volume (p. 39). 

-the specific heat at constant pressure (p. 39 ). 

-a general collision integral (p. 126) (p. 134 ). 


-the space-like element of three volume (p. 20 }). 


-the volume element in momentum space (p. 20). 
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do, -the element of solid angle (p. 93 ). 
-, - mr : 
e a specific combination of Dy Sah Be lope ee 
ne -a coefficient in the expansion of the Boltzmann 
equation for massless particles (p. 105). 


GO! -a coefficient in the expansion of the Boltzmann 


equation for massless particles (p. 145). 


Doha ~the collision term in the Boltzmann equation (p. 
p\ Ee 

en -the electric charge of a particle (p. 19). 

f, -the first order deviation from equilibrium 


expressed aS a power series in a, (IE Le 


g -the determinant of the metric tensor (p. 20 ). 

a the spinewe 1G ac Lol Dae 22"). 

Sop -the metric tensor (p. 17). 

h -Planck’s constant (equals one in this thesis) (p. 
inca 

hy =thesmoment unmet |Uxen pe 5008) \Dee 11. 

hy (v,) -the spectral momentum flux (p. 94 ). 

H -the enthalpy (p. 37 ). 

aa -a standard integral (p. 34 }. 

ae -a standard thermodynamic function of a, and 6, (Pp. 
eh Pe 

i, =the partnclesdritie (po. so0) }. 

1OrD -the spectral particle drift (p. 94). 

ee -a standard thermodynamic function of a, and B, (p. 
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Kee") -a frequency dependent collision integral (p. 
TAO 22 

ee) -a frequency dependent collision integral (p. 
140. dig 

K, (py) -a frequency dependent collision integral (p. 
140 oe 

KA (pe) -a frequency dependent collision integral (p. 
Ad es 

Re. -a frequency dependent collision coefficient (p. 
dee 

= -a frequency averaged collision coefficient (p. 
Wi4n) | 

kp -an arbitrary spatial vector of unit length (p. 
932) 4 

L apna =areo llisioneceefficientsdp. 128): 

Die dco kisi onieocemricvent (ps. 130:)% 

Lae -Arcohiasionsintegqrals(p. el27elm 

Dp gthe restimass of a particle (p., 19). 

Mm, -a generalized mass (p. 20 ). 

Mi, Stewie s Sankt Ucn. 2 om) = 

Me =the zeroth onder sparteor theamass flux (bp. 46) 

ee) -a spectral collision coefficient (p. 143). 

a, -a frequency averaged collision coefficient (p. 
109) 

ny -the mass density (coincides with number density 


for massless particles) (p. 29 ). 


: -the number density (p. 35 ). 
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-the Eckart temperature (p. 99 ). 

-the matter temperature (p. 99 ). 

-a frequency dependent collision term (p. 137). 

-the energy-momentum tensor (p. 26 ). 

-the total energy-momentum tensor of the gas (p. 
26 ), 

-the zeroth order part of the energy-momentum 

tensor (p. 47 ). 

-the spectral energy-momentum tensor (p. 97 ). 

-the thermodynamic functions appearing in the 

transport equations (p. 83 ). 

-the collision term in the Boltzmann equation for 

massless particles due to annihilation and 

creation processes (p. 91 ). 

-the collision term in the Boltzmann equation for 

massive particles due to creation and annihilation 

processes (p27) 92%) . 

-the internal energy (p. 37 ). 

-the unit flow vector or the four velocity of a 

comoving observor (p. 24 ). 

-the unit flow vector for species A (p. 42 ). 

-the thermodynamic functions associated with the 

double momentum flux (p. 54 ). 

-a frequency dependent collision term (p. 137). 

-the double momentum flux (p. 27 ). 


-the zeroth order part of the double momentum flux 


(p. 54). 
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V -an arbitrary parameter; also, a specific volume 
(Cie 37ire 
we “WNeRTOur VElocttyictearpaniicie si e520 aa 


* * 
MED oP |b ,s PEIN Wap -transition probability for 


bainamy cola sionsepe 2261) 


* * 
IN Sea poeta: W ABL =trans itt ons probability for 


creation processes (p. 91 }). 
* * 
Gey Perce phere) Wrap =transi tion propabl lity sor 


annihilation processes (p. 91 ). 


x -the space-time coordinates (p. 18). 

Sey -the natural logarithm of IW G8) eae 

Vp -the natural logarithm of N,/b, (p- Lae 

Oy -the thermal potential (p. 32 ). 

4 -the coefficient of volume expansion (p. 38 }. 

a -the mass times the inverse temperature (p. 32 ). 

Be -the mass times the inverse of the Eckart 
temperature (p. 99 ). 

Bu -the mass times the inverse of the matter 
temperature (p. 99 ). 

B. -the inverse of kT, referred to as the inverse 
temperature (p. 42 ). 

Be -the inverse temperature times the flow vector (p. 

427) 

By -the flow vector times 8, (p. 42 ). 

Y “the adiabatict index! (pig407)). 

on -the difference between u, and a common flow vector - 
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-a fitting and frame change of £, Np 61). 
-the first order part of the mass density (p. 
61)" 
“the first Sorder  partvot the distribution function 
(pee44as 
-the first order part of the mass flux (p. 46 ). 
-the first order part of energy-momentum tensor 
(po. Troe 
-a frame change of th u'(p. 43). 
-the first order part of the double momentum flux 
(Dal cemne 
=the Viracedeltaerunction (bp. 20° ). 
-a fitting change of the thermal Be 
ASG 
-a fitting change Si (pReas @)4 
-the first order part of the energy density (p. 
A Fecha 
-the Bose enhancement or Fermi exclusion function 
FGRECON LPS 1emsmip ! (20) 
-the spatial projection operator of a (Gey Ze ae 
-a fitting and frame change of the distribution 
TUNGET OME Dp earo ley. 
-a variable defined by a contraction of the double 
momentum tensor (p. 53 ). 
-the relativistic enthalpy (p. 38 ). 
-a thermodynamic function defined via a derivative 


of n,with respect to B,\Pe 6F 7. 
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"ee -the Minkowski flat space metric (p. 17 ). 
6 -the volume expansion; an angular variable in 


integrals (p. 24 ). 


8 (x) -the heaviside step function (p. 20). 

O, -the thermal potential per unit mass (p. 36 ). 

K -the isothermal compressibility (p. 38). 

a Tapspectnal coefficient of absorption (p. 145). 

cw -a frequency averaged coefficient of absorption 
(p. 108), 

d -the mean free path (p. 64 ). 

AL -the coefficient of thermal conductivity for 


massless particles (p. 112 ). 


ae -“thesheateconductini tyvmnatr ix) tp. 18e5) . 

Ay -a thermodynamic function (p. 48 ). 

vp -the bulk viscosity for massless particles (p. 
116%) 

‘a -the frequency (p. 93 ). 

oe -the shear viscosity matrix (p. 88 ). 

ve -the shear viscosity for massless particles (p. 
2 

ve =“thesbulk viscosity matrix) lp. 987 a) 

aN -a combination of ¥, and @(p. 23 ). 

3 -ditlmnctionos thesdastribution tunction (pb. 28). 

1 - 3.14159...; the bulk stress of the whole gas (p. 
atthe 

7, -the bulk stress of species A (p. 35 4 

(2) -a projection operator (p. 34). 
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=the energy density (p. 30 }. 


A 
Te, -the spectral energy density (p. 94 ). 
fo) ~the entropy per unit mass, the cross-section (p. 
Ble) 2 
ap 
Oo -the shear tensor (p. 30), 
t -a particle’s world line parameter, the optical 
depthetpue2? je 
or -a variable defined by a contraction of the double 
7 momentum flux (p. 53 ). 
db -an angular variable in integrals (p. 153). 
Ds “a function: solely of the distribution function 
Wee NS) OED 
X -a hyperbolic angle variable in integrals (p. 
53a) 
ae -an arbitrary thermodynamic function of position 
Ue Oa, 
nN sdhmdb br tnaky thermodynamic pelsteisihelal (Vermeulen uh: 
= isi fieq t Ae sien Pe 
ey, a collision coefficient (p ) 
ioe -a collision coefficient (p. 131). 
AB 
XY =a) COl i SionsGochiicient. Woe. 129). 
ABnq 
an -an arbitrary tensor function of position and 
momentum (p. 23). 
wer “the vorticity tensor (p. 30°). 
on -a thermodynamic function (p. 50 ). 
nay -a thermodynamic function associated with the 
all 


solution for the deviations from equilibrium (p. 
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special Symbols 


-the spatial derivative (p. 66 ). 

-the time derivative (p. 66 ). 

-over thermodynmical functions it represent a 
derivative with respect to the inverse 
temperature; over collision coefficients it 


represents a frequency avaerage (p. 66 , p. 109). 
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